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It has been assumed that in the simplest case every 
atom that hits a crystal surface of the same material 
condenses. The corresponding speed of sublimation had 
been calculated theoretically. The formula for the equilib- 
rium pressure is changed in the case of atoms, if the electron 
weight in the gaseous state is different from the solid state. 
The influence of this fact on the speed of sublimation is dis- 
cussed in the first part of the paper and it is shown that for 


true metals the speed of sublimation is probably increased, 
while for non-metals a reflection coefficient exists. For 
molecules the equilibrium pressure is higher than for 
atoms because in sublimation there is a transition from a 
limited oscillation of the axes to free rotation. It is made 
plausible that this results in an increased speed of sublima- 
tion but a finite reflection coefficient might also occur. 





TATISTICAL mechanics determines uniquely 
the vapor pressure or the equilibrium density 
of the saturated vapor over a crystal. On the 
other hand, equilibrium is present if the number 
of evaporating molecules is equal to that of the 
condensing. This relation is discussed in the 
following for several cases. It has a close con- 
nection with two recent papers that correlate the 
reaction velocity with the type of bond affected 
by the reaction.! 


I. ATOMS. 


Given a crystal made up of atoms thut have 
no resultant moment 7; let Qo be the extrapolated 
heat of sublimation at zero temperature. Let 
furthermore vy be the (geometrical) average of 
the frequencies of vibration present in the 
crystal, then if the temperature is high enough 
to have classical specific heat 3R, the number of 
vapor atoms per cc Np is given by? 


No=(24M/RT)'v* exp (—Qo/RT)=Ko. (1) 


'O. K. Rice and H. Gershinowitz, J. Chem. Phys. 2, 
853 (1934); H. Eyring, J. Chem. Phys. 3, 107 (1935). 

*K. F. Herzfeld, Kinetische Theorie der Warme, Braun- 
schweig, 1925, p. 223. 


If the temperature is so low that the specific 
heat of the solid is zero, 


No=(2eMRT)\Lh)- exp (—Q:/RT)=Kv' (2) 


(where L is the number of molecules per mole 
and Q; is Q) minus the zero point energy). In 
equilibrium the number of molecules condensing 
is? 

(RT/2xM)'No (3) 
if one assumes that all atoms hitting the surface 
actually condense.‘ If we designate the number 
of atoms evaporating under this condition per 
second and cm? with kp and abbreviate 


w=(RT/2rM)}, (3’) 


the x component of the average thermal velocity 
in the gas, there is in equilibrium 


3]. Langmuir, Phys. Rev. 8, 149 (1916). 

‘ This is in agreement with experiments on Hg, Cd, Zn 
by M. Knudsen, Ann. d. Physik 50, 472 (1916); and 
R. W. Wood, Phil. Mag. 30, 300 (1915); 32, 364 (1916). 

5 The formula for k) which results from (4) has been 
deduced from kinetical considerations directly on p. 232 of 
book cited, footnote 2; see also R. Becker, Zeits. f. Physik 
28, 256 (1924). 
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ko=wKo, (4) 


ko=wNo or 


or, at low temperatures 


ko’ =wKy’. (4’) 
Assume now we have an atom that has in the 
gaseous state a moment j, therefore an electron 


weight 
£v - 2j+ 1 ’ 


while the crystal is not distinctly paramagnetic 
(this includes metals) so that the electron weight 
in the crystal is 1. 

Then we have to replace (1), (2) by 


N=Kyog, or N=Ko'g». (5) 


If we call now the speed of evaporation k or k’, 
we find 
(6) 


where 1—a is the reflection coefficient, i.e., the 
fraction of the atoms that hit the surface but 
do not stick. Defining now a number 6, k=), 
we find 

(7) 


k=wNa, 


B/a=». 
The discussion is mainly concerned with (7). 


(a) Emission of free electrons from a metal 


Here g,=2. The direct calculations of Nord- 
heim® show that a=1, B= 2. The physical reason 
for this is that an electron entering the metal has 
such a high energy that it gets into an empty 
part of the energy band of the metal and can 
therefore enter, independent of the direction of 
its spin. On the other hand, there are, because 
of the spin, twice as many electrons of a given 
state in the metal than without spin, and there- 
fore B=2. 


(b) Metal with one valence electron in Si/2 


Here again g,=2, because of the two possible 
spin orientations, while the experimental fact that 
there is no strong, temperature dependent para- 
magnetism shows that we need not consider 
different orientations in the crystal. Two cases 
must now be distinguished. 


(b’) 
If the solid is a true metal and can therefore 
be considered according to Wigner, Seitz and 


6L. Nordheim, Zeits. f. Physik 46, 833 (1928); Proc. 
Roy. Soc. Al21, 626 (1928). 
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Slater’ as a lattice of positive ions surrounded 
by a sea of Fermi-electrons, one should expect 
a=1. Each condensing atom would create two 
new neighboring states of opposite spin and 
therefore the electron of the condensing atom, 
in dissociating off, would find a ready free level. 
The ion, of course, having j7=0, would always fit. 
This result is in agreement with Knudsen’s 
experiments, who found that every atom of Ag 
and Cu condensed.* One must, however, account 
then for B= 2 as a consequence. 

The electron of a condensing atom will in 
general get into the free part of the energy band 
and drop down later. Therefore, from the 
principle of detailed balancing, one would expect 
that an ion ready to evaporate would in general 
pick up one of the few electrons in the unfilled 
part of the band. There will be accordingly 
certain electron states of energy Q)’/L in 
which the electron ‘‘belongs” to a certain ion. 
The fraction of the ions that will have electrons 
in these states will be 


2e~2o'/RT, (8) 


the factor 2 coming from the two possible spin 
directions. The number of ions having such 
electrons that evaporate per second is then given 
by the formula (4) in which Qo)’=Q)—Q)’ is 
inserted, multiplied with (8), which gives p=2. 
If the one valence electron is in a state different 
from S in the vapor, while the conditions in the 
solid are as described here, the above considera- 
tions hold with 


go=2j+l1, B=2j+1. 
(b’’) 

If, on the other hand,° it were possible to 
condense these atoms into a monatomic, non- 
metallic, diamagnetic crystal, in which, therefore, 
the atoms keep their valence electrons and have 
orientations of the moment fixed in the lattice, 
g, would still be 2 (or 27+1 in general), g,, the 
weight in the solid, 1. But now the considerations 
that led to the deduction of the evaporation 


7J. C. Slater, Rev. Mod. Phys. 6, 209 (1934). 

8 M. Knudsen, Ann. d. Physik 50, 472 (1916). : 

* A distinction analogous to that between (b’) and (b”) 
is also made by J. H. Van Vleck, The Theory of Electric 
— Magnetic Susceptibilities, Oxford, 1932, p. 283, cases 
e, f. 
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formula (4) are valid without change, therefore 
§=1 and @ must be 3 (or 1/(2j+1) in general). 
This can easily be understood. If an atom comes 
down from the gas, only half (or in general the 
fraction 1/(27+1)) of the places on the crystal 
on which it can alight have as permitted direction 
of moment the one which the incoming atom 
has; this results in a reflection coefficient of the 
magnitude given above. Formally one can com- 
pare (b’) and (b’’) so: A condensing atom with 
one valence electron creates two new states in 
the electron term-scheme of the crystal. But 
in the case (b’) the two levels corresponding to 
the two spin directions are very close together, 
in the case (b’”’) rather far (i.e., much farther 
than kT) apart. 


(b’’’) 


If the crystal should be strongly paramagnetic, 
so that the 7 in both states would be the same 
and the components in the direction of a mag- 
netic field the same, g,=g;,, a=@6=1. 


II. MOLECULES 


We assume that the molecule in the crystal 
has, besides the vibrations of its center, rotatory 
oscillations of its axes (or axis for a linear 
molecule) around an equilibrium position.!° We 
assume the temperature to be so high that these 
oscillations have their classical specific heat 
(3R, or 2R for a linear molecule), and that the 
specific heat of rotation in the gas has the 
classical specific heat ((3/2)R, or R for a linear 
molecule). Finally we assume that there is 
either no difference between the internal vibra- 
tions of the molecule in the gas and in the solid 
(the »’s remaining unchanged) or that in both 
states they do not contribute to the specific heat. 

The vapor pressure formula is then 


N= KoK,. (9) 


Here K, is the ratio of the phase integrals of the 
rotational coordinates in the gaseous and the 
solid state. In classical calculations the phase 
integral can be split into a factor depending 


"D. H. Andrews, Proc. Amst. Acad. 29, 744 (1926); 
see also A. Eucken, Jahrb. f. Rad. 16, 364 (1920); A. 
Eucken and E. Karwat, Zeits. f. physik. Chemie 112, 479 
(1924); A. Eucken and F. Fried, Zeits. f. Physik 32, 151 
(1925); F. Simon, ibid. 31, 225 (1925). 
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only on kinetic energy and momenta and one 
depending only on potential energy and co- 
ordinates. The former is independent of the 
forces and therefore the same for the gas and the 
crystal. Accordingly it drops out of K,, which 
can then be written 


K,=%,°?,"", 


(10) 


where ® is only the integral over the coordinate 
space. 
But we have always 


o>. 


These quantities are, namely, the average steric 
angles covered by the motion and these are 
always much smaller for the crystal than for the 
unhampered rotation of the gas. 

For example, for a diatomic molecule with 
unequal atoms (symmetry number 1), ®, is 


simply the full spherical angle 47, while 
&,= fdef exp (—ILw*d*?/2RT) sin ddd, (11) 


where I is the moment of inertia, w the average 
frequency of oscillation. This gives (putting 3 


for sin #) 
©,=2rRT/IwL, 


or, e.g., with w=3.10", J=10-**, T= 300. 
},=47-0.02 


or K,~50. 
We have then in general 


B/a=,?,"! (12) 
and again the question arises: 

(13’) 
(13’”) 


is B=,?,—', i.e., very great, and a=1 


or B=1, a=,%,"', i.e., very small 


or does something between these extremes 
happen? Experimental results seem not available. 

We are going to show that (13’) and (13’’) 
follow clearly from two different assumptions. 
The assumption leading to (13’) seems more 
probable. 

The question is the following. The heat of 
sublimation per mole, extrapolated to T=0, is Qo. 
Consider a molecule in its heat motion. Assume 
that in a given moment it has a potential energy 
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(calculated from the equilibrium position) of e, 
which is made of two parts, e, which is due to 
the deviations of the center of gravity of the 
molecule from its equilibrium position, and e,, 
which is due to a rotation of the axes out of the 
equilibrium orientation. If now the kinetic 
energy necessary to tear the molecule loose is 
given by 


e’ =Qo/L—€=Qo/L—e:—e,, 


that is, if the potential energy stored up by the 
oscillation of the axes is available for sublimation, 
(13’) holds. If, on the other hand, it is not 
available, that is, if the energy necessary to 
sublimate this molecule is 


e’ =Q)/L—e: 


(14’) 


(14”) 


(13’’) is true. 

To prove that, start with the latter, simpler 
case. Then the deduction given for k= kp can be 
taken over without change, B=1. a=,%,— can 
then be understood so that only such molecules 
can be condensed which come down under such 
an angle of the axes with the equilibrium position 
as it occurs at the temperature J in the crystal. 
(This is, of course, to be understood in the 
proper way, i.e., considering the distribution 
law.) The fraction of these properly oriented 
molecules is obviously given by the value of a 
mentioned above. 

If, on the other hand, (14’) is true, we have tc 
modify the deduction leading to (4) in the 
following manner: 

Call g in general the coordinates that char- 
acterize the orientation of the axes, dg the 
corresponding element in the phase space. Then 
the fraction of the molecules that have deviations 
from the position of equilibrium between x and 
x+dx--- and a direction of the axes between 
gq and g+dq is 
v3(2xM/RT)3e—"«"®"dxdydz 

Xb, le-LelkTdg, (15) 


The number of molecules hitting the surface 
(y, 2 plane) so that mz?/2=Q)/L—e is 


n'(M/2eRT)'f* | exp (—Ma*/2RT) de 


=n! (RT/24M)ie-(%-L9/2T, (16) 
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where 7’ is the number of the particular molecules 
per cc. Inserting (15) into (16) and integrating 
over dxdydz and dq, one finds as number of 
evaporating molecules 


v3(2nM/RT)! exp —(Qo/RT)(RT/2xM)'@,-'fdg 

=kyb,'fdg. (17) 
Sdq=®,, (18) 
B=0,0,-1. 


But 


and therefore 


In (18) it is understood, that the integration is 
to be taken according to classical physics, so 
that for example for a diatomic molecule with 
equal atoms, symmetry number 2, 


Sdq=227=47/2. 


In that case, a=1. To understand why under 
these assumptions even badly oriented molecules 
condense, one has to reverse the reasoning used 
for (15’). The large forces involved in the large 
heat of condensation Q/Z are sufficient to 
force the molecule into the lattice in whatever 
orientation it arrives, providing the necessary 
potential energy e,. 

Under both assumptions (14’) and (14’’), the 
kinetic energy of the oscillation of the axes is 
not involved in the process of sublimation. That 
means that each molecule keeps the angular 
momentum it had just before sublimation also 
in the gaseous state and each molecule keeps the 
angular momentum in the moment of condensa- 
tion it had just before. Of course, angular 
momentum begins to change immediately after 
condensation, because of the forces now acting 
on the molecular orientation. 

If the temperature rises so high that the 


molecules rotate almost freely in the crystal, 
®,~®,, B=a=1. 


III. REFINEMENT OF THE THEORY 


In the preceding deductions, one assumption 
has been made that might be wrong. The 
reasoning leading to (4) said: The number of 
atoms that have a deviation from the equilibrium 


position between x and x+dx--- and evaporate 


per second is 


n'(RT/2nM)} exp (—Qo+Le/RT), (19) 
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where €, is the potential energy connected with 
that deviation and n’ the number of atoms per 
cc having this deviation. One then inserted n’ 
from the Maxwell-Boltzmann distribution law 
and integrated over dxdydz. But that is correct 
only if in the permitted range Le,<@Qp. Other- 
wise, One can only integrate to such values of 
x, y, in (19) that 


Le = Qo. 


Call Sok*—%dxdydz=1/n""’. 


Then, in the simplest case, the number of 
evaporating atoms is 


nko/n" =a ko, 


where ” is, as before, the number of atoms per 
cc, and speed of sublimation and speed of 
condensation have both to be multiplied by 
a’<1 in all the cases treated in I. The same is 
true for molecules and assumption (14’’). Possi- 
bly Volmer’s" experiments showing a reflection 
coefficient of Hg on mercury crystals of 10 
percent can be explained so. 


1M. Volmer and I, Estermann, Zeits. f. Physik 7, 1 
(1921). 
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Similarly for molecules, we have to change the 
discussion of assumption (14’) so that we inte- 
grate only to such values of the coordinates in 
the crystal that satisfy 


Qo>Ll(ete,). 


The speed of sublimation k8 has again to be 
multiplied by a factor a’ smaller than 1, 


a! =n®,—f 5% xdydadq =’, 


and the same factor enters the formula for the 
speed of condensation, giving rise to a reflection 
coefficient of 1—a’. 

The physical reason for the existence of this 
reflection coefficient is, of course, that the 
attractive forces that are responsible for the heat 
of sublimination can force the incoming molecule 
only into the fraction a’ of the total coordinate 
space. 

We have then an intermediate case between 


(13’) and (13’) where 
k=k@', B=8'O-1, of =8'8.—, 


I acknowledge gratefully the clarification of 
the problem due to discussions with Dr. M. 
Goeppert-Mayer. 
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Non-Silicates with Cristobalite-Like Structure 


Tom. F. W. BArTH, Geophysical Laboratory, Carnegie Institution of Washington 
(Received March 30, 1935) 


The cristobalite type of structure which has been previously observed in certain silicates: 
carnegieite, NaAlSiO,, and NagCaSiO, is now shown to occur in K,0.Al,0;, K2O.Fe.O;, and 


Na.O. Al.Os. 


NEW cubic compound, K,O.Al,O;, has 
recently been described by Brownmiller,! 
and measurements on an x-ray powder photo- 
graph of it have been published.? This powder 


‘L. T. Brownmiller, A Study of the System Lime-Potash- 
Alumina, Am. J. Sci. 29, 260 (1935). 

*The reflections observed by Brownmiller correspond 
to cubic symmetry (length of the edge of the unit cube, 
®=7.69A), except for three weak lines, corresponding to a 
Planar spacing of 2.395, 2.000 and 1.488. He states, 
however, that K:0.Al,0; is extremely hygroscopi¢ so 
there were considerable difficulties in taking the photo- 
graph; it is probable, therefore, that these weak reflections 


— - some alteration product on the surface of the 
rystals, 


photograph is very similar to that of KO. Fe:Os;, 
and they both resemble the x-ray patterns of 
certain silicates structurally isomorphous with 
cristobalite (see Table I). These similarities 
indicate that the structure of K,O.Al.O; and of 
K;,0.Fe,O3; may in some way or other be related 
to that of cristobalite. 

The well-known feature of the cristobalite 
structure is the existence of endless silicon-oxygen 
chains with large, open spaces in between. It 
seemed probable, therefore, that in K,O.Al,O; 
there would be corresponding chains of alumi- 





T. F. W. BARTH 








Fic. 1. The unit of the high-cristobalite structure. The 
large transparent balls represent oxygen ions, and the 
small black ones silicon ions, bringing out their relative 
size. The positions of the ions are very slightly displaced 
to give the structure a more regular arrangement than 
actually exists. 


num-oxygen, and in K,0.Fe.O; chains of iron- 
oxygen, with the potassium atoms situated in 
the interstices.* 


TABLE I. Intensity of powder diffraction lines of a-carnegieite, 
Na,CaS10,, K2FeO,4 and of K2Al04. 
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* Observed intensities as given by Brownmiller. 


3 For a description of the cristobalite structure, and of 
the structure of carnegieite and Na,CaSiO,, which have 








Fic. 2. The unit of the structure of KsAl.O, (and of 
KiFeiO.), The large transparent balls represent oxygen 
ions, and the small black ones aluminum ions, bringing 
out their relative size. In this way aluminum-oxygen 
chains are formed that are strictly analogous to the 
silicon-oxygen chains of cristobalite. The centers of the 
potassium atoms are indicated by crosses; the atoms 
themselves are so large as to fill up the interstices, touching 
12 adjacent oxygens. The positions of the ions are very 
slightly displaced to give the structure a more regular 
arrangement than actually exists. 


The constituent atoms in cristobalite (and in 
the silicates isomorphous with it) show a close 
approach to a face-centered arrangement with 
silicons close to 4b, 4d; and oxygens close to 16a.‘ 
In K,O.Al,O; one would expect, therefore, Al to 
be in 4b, 4d; and the oxygens to be in 16a, or 
very nearly in these positions. The only possible 
place for the K-atoms would then be 4c, 4¢. 
The same holds, of course, for KO. Fe:O3 with 
interchange of Fe and Al. 

We can test this assumption by calculating 
the intensities which, if this represents the true 
structure of K,0.AI,03; and K,O.Fe.O3, must be 
proportional to the observed ones. 

Table I shows that there is a fair agreement 
between calculated and observed intensities, and 
it can be concluded, therefore, that the structure 


the‘same type of structure, see Barth and Posnjak, Silicale 
Structures of the Cristobalite Type: III, Structural Relation- 
ship of High-Cristobalite, a-Carnegieite, and Naz ‘a Si04, 
Zeits. f. Krist. 81, 376 (1932). 

4 Wyckoff’s notation. 
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of K,0.Al,0O3 and of K,O. FeO; is of the cristo- 
balite type. 

A few years ago it was shown that the structure 
of cristobalite was not exactly face-centered but 
actually built upon a simple lattice.5 Also certain 
silicates isomorphous with cristobalite exhibit 
the same small, but characteristic deviation from 
a face-centered lattice.® 

An inspection of Table I shows that KO. Al.O3; 
and K,O. Fe:O; also seem to deviate slightly from 
the face-centered arrangement. The discrepancies 
between observed and calculated intensities are 
small, but significant ; thus 400 is always stronger 
than demanded by the theory, which is again 
exactly the same phenomenon previously ob- 
served in cristobalite and its congeners. This 
means that the structural arrangement we have 
been considering so far is not accurate, it 
represents only the approximate or idealized 
structure, whereas the true structure is slightly 
different and actually isomorphous with the 
space group 7‘. (The arguments for 7“ are ex- 
actly the same as in the case of cristobalite.) 

From this it follows that the true structure of 
KO. Al,O3 is conditioned by no less than 8 inde- 
pendent parameter values: 4K+4K, 4Al+4Al, 
and 40, in position 4f, each group of four with 
one parameter, and the remaining 120, in the 
general positions of 7* with three parameters. 
These 8 parameters cannot be determined accu- 
rately from powder photographs. But in the case 
of K,O. FeO; a trial calculation has been made 
which shows that the discrepancies between 
calculated and observed intensities are appreci- 
ably diminished by moving the oxygen atoms 
slightly out of the idealized positions, following 
thereby exactly the same procedure previously 
described in the case of carnegieite and NaeCa- 
Si0,. 

Figs. 1 and 2 show very plainly the relation 
between the structure of cristobalite and the 
Structure of KO. Al,O; (or K.O. Fe.Qs). Al (and 
Fe) has coordination number 4, K has 12. 

The interatomic distances in the idealized 
structures are as follows: 

*T. F. W. Barth, Am. J. Sci. 23, 350 (1932). 


‘T. F. W. Barth and E. Posnjak, Zeits. f. Krist. 81, 
135, 371 (1932). 


CRISTOBALITE STRUCTURE 


For KO. Al,O;: O —O=2.72A, 
Al —O=1.66A, 
K —O=3.19A; 
O —O=2.82A, 
Fe—O=1.73A, 
K —O=3.32A. 


For K.O f FeO; ° 


It is notewerthy that the distances Al—O and 
Fe—O are smaller than usual even for coordina- 
tion number 4,’ a phenomenon which also has 
its strict analogy in the silicate structures of the 
cristobalite type. 

As has long been known, mixed crystals of 
a-carnegieite and cristobalite exist. But it has 
now been observed that a-carnegieite is able to 
take up also Na,O.Al.O; in solid solution; 
indeed, it is possible that a continuous series of 
solid solutions extends all the way between 
carnegieite and Na,O.Al.0;.° Consequently 
Na.O. Al,O; must be isomorphous with carnegieite 
and the structure type of at least one of its 
modifications (potential or real) must therefore 
be that of K20.Al,03;. This becomes so much the 
more probable since powder photographs of 
Na,O.Al,O; at ordinary temperature exhibit an 
interesting similarity to those of low-carnegieite. 
To be sure, the low-carnegieite pattern merges 
gradually into that of low-Na,O.Al.O; through 
intermediate mix-crystals with successively de- 
creasing amounts of silica.'® 

Thus we see that the cristobalite type of 
structure is not restricted to silicates only, but 
that at least two aluminates and one ferrite 
crystallize in the same structure. 


Note added in proof: In a paper by S. Hilpert and A. 
Lindner measurements of x-ray powder photographs of 
K,0.Fe0;, Rb:O. Fe.0;, and PbO. Fe,0; have been 
published (Zeits. f. physik, Chemie B22, 400 (1933)). From 
these measurements it can be deduced that Rb,O. FeO; 
and PbO.Fe,0; are isomorphous with K,0.Fe.,0;, and 
consequently that these two additional ferrites also 
crystallize in the cristobalite type of structure. 


7 For discussion and values of interatomic distances in 
relation to coordination number see W. H. Zachariasen, 
Zeits. f. Krist. 80, 137 (1931). 

8 Cf. Barth and Posnjak, Silicate Structures of the Cristo- 
balite Type: ITI, Zeits. f. Krist. 81, 385 (1932). 

® Unpublished experiments by J. F. Schairer of this 


‘Laboratory. 


-10 These mixtures were prepared by Dr. J. F. Schairer. 
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The Moles of Vibration of Butane and Pentane. “Free Rotation’’ about Carbon- 
Carbon Bonds and a New Type of Stereoisomerism! 


Louts S. KAssEt,? U. S. Bureau of Mines 
(Received March 16, 1935) 


The Hamiltonian function for a chain of four carbon 
atoms with the simple valence force system potential 
energy is transformed by successive changes of variables 
to such a form that the wave equation is nearly separable 
into five harmonic oscillator equations and a sixth equation 
governing the internal rotation. This rotation equation 
has been solved numerically for twelve different rotational 
levels of the nonvibrating molecule, some of which actually 
correspond to rotations, others to a torsional oscillation of 
restricted amplitude about either the cis or the trans 
position. The remaining terms in the wave equation are 
treated as a perturbation; the first-order perturbation 


energy is less than 1 cm™ for all moderately low vibrational 
states. There are thus two isomeric forms of butane with 
plane carbon skeletons, and a third form in which the two 
ethyl groups are rotating about the central carbon-carbon 
bond. The vibration frequencies for the three possible 
plane forms of a five-carbon chain have been calculated. 
When the frequencies of C2, C3, Cs and C; are calculated 
with the same two force constants throughout, and by 
taking account of the two plane forms of C, and the three 
of C;, results are obtained in fair statistical agreement 
with observed Raman and infrared lines. 





INTRODUCTION 


ERY little work has been done on the 

motions of molecules which are capable of 
internal rotations. Nielsen’ has studied a general- 
ization of the symmetrical top, one with an 
extra degree of torsional freedom about the 
symmetry axis, and LaCoste‘ has treated the 
spherical top with four symmetric tetrahedrally 
arranged gyroscopes. These investigations of 
models of ethane and tetramethylmethane, re- 
spectively, both deal with the interaction be- 
tween external and internal rotation for non- 
vibrating molecules. The present paper is an 
attempt to consider the interaction of the in- 
ternal rotation with vibration. This interaction 
will be relatively small for the foregoing mole- 
cules; for saturated hydrocarbons with a chain 
of four or more carbon atoms, the internal 
rotation changes the shape of the carbon skele- 
ton, and a relatively large interaction results. 
Because of the considerable difference in mass 
between carbon and hydrogen atoms, there is 
little interaction between C—H and C—C vibra- 
tions, and useful results probably can be ob- 
tained by treating the carbon skeleton alone. 
In addition to neglect of the hydrogen atoms, and 


1 Published by permission of the Director, U. S. Burea 
of Mines. (Not subject to copyright.) . 

2 Associate physical chemist, U. S. Bureau of Mines, 
Pittsburgh Experimental Station, Pittsburgh, Pa. 

3 Nielsen, Phys. Rev. 40, 445 (1932). 

4 LaCoste, Phys. Rev. 46, 718 (1934). 


treatment only of states in which the total 
angular momentum of the molecule is zero, it is 
necessary to use the simplest possible valence- 
force potential function 


V=(k/2)D Ad? +(b/2)Xd*A0*, (1) 


where d; and 6; are the bond lengths and bond 
angles; d, the equilibrium bond length, has been 
introduced in the second term to make k and } 
dimensionally homogeneous. With these three 
simplifications, it becomes possible to give a 
fairly complete solution of the problem. 

For completeness, we give the results for C: 
and C; based on the foregoing assumptions. In 
terms of \ defined by 


vy=(1/27)d1/? 
they are 


Co: \=2k/m, 
Cs: X=(7/3)(k/m), 
d= (5/3)(k/m) + (16/15)(b/m) 
+ (288/125)(b?/km)+--- 
A= (18/5)(b/m) 
— (288/125)(b?/km)+-:-: 
the last two values being roots of 


d?—(5k+14b)(A/3m) + 6bk/m? =0. 
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HAMILTONIAN FOR THE C, CHAIN 


The coordinates in which the potential energy 
is written are shown in Fig. 1. It is desirable to 
make the substitutions 


ditd3=de,, de=de2,, di—d3=de3, (8) 
6:+62=m1, 01—02= 2. (9) 


The kinetic energy is most easily obtained from 
auxiliary Cartesian coordinates referred to the 
center of gravity as origin, thus excluding trans- 
lational energy. It is not feasible to exclude 
rotational energy of the whole molecule as can 
be done in simpler problems by using an auxiliary Fic. 1. 

system which rotates with the molecule. The 

method used here should be generally useful for first calculate the kinetic energy of vibration 
polyatomic molecules of low symmetry. We plus rotation 
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4 
T’=(m/2)>- (é2+92+ 2,2) = (md?/2)[(9 —4y?) €12/18 + €1€2/3 +2 ?y7€191/9+ 8y5€19/94+ eo” 
i=1 
+ 23!€oqi + (6 —-y?) 917/12 —29?y5919/9+8(1 —y?) G?/9 ] 
+ (md?/2)((9+8y7) €3?/36 — 2"/*y7e32/9-+ (4+-y*) j2/36], (10) 


where y=cosy, 6=sin ¢. (11) 


This expression is not exact, having been obtained by putting d:\=d.=d3=d and 0,= 02=cos™! (—1/3) 


wherever functions of these variables occur as coefficients in %;, ¥; and 2;, in accordance with the 
usual custom in treating small vibrations. No similar substitution is possible for ¢, and as a result 
the quadratic form for J’ does not have constant coefficients. The error introduced at this stage 
could be reduced by the perturbation-theory methods recently used by Bonner® for the water mole- 
cule, but this refinement would be wasted in the present crude treatment. 

The next step is to determine the external rotational energy of the whole molecule, which must be 
subtracted from (10). The moments of inertia are 


A=m>(y2+2;*) =(md?/9)(16—8y?);  B=my(x?+2,7) = (md*/9) (33 — 167?) ; 
C=m>d(x2+y.2) = (ma?/9)(17+8y?) 


(12) 


and there is one nonvanishing product of inertia 
E=(10/9)2"/?md*6. 
The angular momenta about the Cartesian axes are 
M,=(md?/9)[4ybe3—2"!?y5H2], My=(md?/9)[3-2"76(e1—2€2) — (21/2) dq +5-2"? ye], 
M.=(md?/9)[ —4-2"/*yé3+ (13/2) 72]. 
It is shown in the Appendix that the rotational energy is 
Trot = 3 M,?/B+(CM2+2EM.M.+4M,’)/(AC—E?’) ]. 


When this rotational energy is subtracted, the kinetic energy of vibration is found to be 


Teh 


* Bonner, Phys. Rev. 46, 458 (1934). 
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T =[md?/(297 —144y?) ][(261/4 —60y?+ 167‘) €:?+ (171/2 — 60y”) €1€2 +2"/2(63/2 —15y?—8y*)éim 
+-5(72 — 6472) 6+ (225 /2 — 3672) &?-+ 2236 +152) eo + 120 y5é29-+ (153/8+15y2+274) Hr 
+212 5(72+ 16y*) m19+ (132 —32y?— 6474) o? ]+Lm/(72+192y?—64y*) JL(9+ 2477) és? 

+12-2"?y%€sH2+972"]. (16) 

The potential energy in these coordinates is 


V=[(kde:2/4+hAe?/2+bAn12/4) + (kAes?/4+bAn2?/4) ]d?. (17) 


If y~ were not a function of y, one would now find by straightforward calculation in the usual way 


that the frequencies were (1/27)\"/? where the A were roots of 


d h?—[2k+ (2/3)(3+8y7)b ](A/m) + (4/9) (9 +247? —8y*) (bk/m?) =0 (18) 
an 


h4—[4k+ (22 —16y7)b/3 ]d3/m+[34k?/9+ (8/9) (27 — 16y?—4y*) bk ]A?/m? 
— (4/9) (33 —8y?—16y')bk?h/m>=0. (19) 


These roots may be expressed conveniently as power series in b/k, which is about 0.1; they then 


become 
hi =2k/m+ (16/9) y4b/m+ (16/9) 741+ (8/3)? — (8/9) y* ]b?/km+-:--, 





_ he = 2[1+ (8/3)y?— (8/9) 4 ]b/m— (16/9) y4[1 + (8/3) 7? — (8/9) 74 ]b?/km + - - « ” 
hg = (2+2"2/3)k/m+(4/51)[(22 —28y2+ 1274) —2"/2(15 —16y2+2y4) ]b/m+---, 
Na=(2—2"?/3)k/m-+ (4/51) (22 — 28y?+ 124) +2"/7(15 —16y?+274) ]b/m+---, (1) 
As = (2/17) (33 —8y?—16y4)b/m— (16/4913) (33 —8y? — 16y*) (81 — 1007?+38y')b?/km+---, 
Ae=0. 

The linear transformation which would reduce (16) and (17) to 

T= (ma*/2)) 5 iia mel, (22) 
i=! 9 33-167’ 
V= (md? /2)¥- hu; (23) 

" i=1 

Aes= — 2+ (4/9)2'%42(9-+24y2— By!) 12(b/k) 2, 9 

Ane= (4/3)2"/?y2u1+ (2/3) (9+24y?—8y4) "U0, (24b) 

Aer = (2+21?/3) "ug (2—2"?/3) "ut (4/17) (7 —6y?)[(2/17) (33 — By? — 164) ]"(b/k)us, (24) 

Aeo=(1+2"/?/6)"/?us+ (1 —2"/?/6) "4+ (4/17)(4—?) (2/17) (33 —8y? — 16y*) ]!7(b/k) us, (24d) 


Ani=[(4—25) +297? ][(1/17)(1— 22/6) us 

+[— (44252) +2522] (1/17) (1242/6) }¥a4-+[(4/17)(33—8y?— 1674) "2s, (246) 
g=y6L(1/17) (1 —21/2/6) ]"/?u3+ yd (1/17) (1+21/2/6) |x, 

+276(9+2y7)[(17/2) (33 —8y? —16y*) -Y?us+u6, 
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where every coefficient is merely the leading term of a power series in b/k. The substitution (24) 
does still lead to (22) and (23), even though 7 is a function of the old variable ¢, but of course the old 
variables have not been completely eliminated. We note, however, that uw is the only one of the new 
variables which is not associated with a small amplitude vibration. We can therefore express y 
explicitly by means of power series in the u; 


y =cos uf 1+ > a;juujt+:- [sin uf > buit-:- +} (25) 
i, j=3 i=3 


The coefficients a;;, b;, etc., could be evaluated by substitution of (25) into (24); when (25) is sub- 
stituted into (22) and (23), however, one sees that the 5; occur only in terms of at least the third 
degree, and the a;; in terms of at least the fourth degree. Since third-degree terms have been neglected 
from the beginning in both 7 and V, we continue to neglect them here and adopt the new definition 


YY =COS Ug. (26) 


In terms of the momenta, (22) becomes 





9 33-167? 
T = (1/2md? x p?+- ‘| (27) 


De 
i=: 8 33 —8y?— 1674 


We have now obtained the Hamiltonian in a form for which the variables are nearly separated. 


WAVE-MECHANICAL TREATMENT OF THE Cy CHAIN 
Writing the wave equation by the method of Condon and Morse® and making the substitutions 
ug=(h/2r)"!2(mkd*)-"2n,, i=1,2,3,4,53 we=v6e, (1/27)(k/m)"?=%, Ai=(Rk/m)p;i (28) 
ay ay 


5 2W 
we obtain > —+A—+B—+ -£ pv? Jy=0, (29) 


1 Ov,” OV_62 Os 


where 
A =(9h/322?md?) (33 —16y?) /(33 —8y?— 167%), (30) 
0 
B= (27h/42?md?)76(33 —132y?+32y*)/(33 —8y?— 167‘). 
We shall solve this equation by a perturbation method. Let 


0°*y;/d02°+(2Wi/hvo— piv.2)Wi=0, i=1, 2, 3, 4, 3, (31) 
AdWe/d062-+Bdve/dv6-+ (2/hn)(We—X Wivo=0. (32) 
i=1 


‘ Since the p; in (31) are functions of y and therefore of vs, the W; are also functions of vs, and thus 
play the role of a potential energy in (32). It is then easily seen that 


Yo= Visa e (33) 
5 OW O7Yo Ovo | 2Ws 5 b 


satisfies 


> —+A—+B—4+| —-Y pa?+V! 


i=1 Ov,” OV¢" dvg Lhy = 


5 10%; 1 des S41 OW 1 OY; 8 1 OW 
—Vis Ay ——+2 ene ane a )=3=—luh7 
Rae ee - i= ly; OV_62 Ve OVs i= yp; OV, i=j+1 j= ly; Vs V; Os i=1 Vi Vs 


*Condon and Morse, Quantum Mechanics, McGraw-Hill, 1929, p. 15. 


(34) 


where 
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Then if Yo is taken as the zeroth-order solution of (29), the perturbation potential is evidently V’, 
and the first-order perturbation energy is 


awiyin={ ff ff [vevar. (36) 


Now since f y7dv;=1 (37) 
we find by differentiation with respect to v¢ 
© Oy; dpi\? 7? fOVi? n?+ni+1 /dpi\? 
—dv;,= -( ) f (~) dv;= — ~~ ( ) (¢=1, 2,3, 4,5), (38) 
7 ow Ov,6" dv¢5 —o Op: 32p;? dvs 


where the last result makes use of the fact that the y; are normalized solutions of (31) and therefore 
are 
























Y= p,"8a-14(2im; !)-V2 exp (—2-!p;"/20,2) 17 n,( pi"! 0;) (i=1---5) (39) 





and uses also various well-known properties of the Hermitian polynomials. Most of the terms in (36) 
vanish and the first-order perturbation energy reduces to 


2W! 2s Aye? 5 dp; 2 
= -{ 32 2 (ni? +ni+1)pi- ( ) dv¢. (40) 
0 


hv dvs 











The zeroth-order energy, of course, is 


2W?/hy=2We/hv. (41) 


















This is as far as the analytic solution can be _ tions of (32) subject to the boundary condition 
carried ; (32) is not amenable to anything except 
: . : = , 42 
numerical solution. Fig. 2 and Table I show the velo) delet 2m) (42) 
dependence of the p,/? upon g, calculated 
numerically from (18) and (19) for the case 
5 


b/k=0.1, and also the composite (1/2)>°p;!/? 
t=1 


which is the potential energy term for (32) 
when all five vibrational quantum numbers are 
zero. This curve has a period of 7, and is sym- 
metric about g=0, 7/2, a---+. The differential 
terms in (32) likewise are seen to be symmetric 
about g=0, 7/2, r---+. The internal configura- 
tion of the atoms has a period of z, but a change 
of 27 is required to bring the molecule back to 
its original position. We therefore require solu- 














RELATIVE VALUES OF Vp 


TABLE I. Values of p;'!* for Cy with b/k=0.1. 














pii/2 





p2i/2 p3i/2 pal/2 psi/2 (1/2) 


1.49550 0.70484 1.58475 1.26027 0.31667 2.68101 
1$° 1.48503 -70165 1.58275 1.26572 -35853 2.69684 
30° 1.45968 -68508 1.57800 1.28387 -44148 2.72405 
45° 1.43385 -64864 1.57347 1.31594 -50908 2.74049 
60° 1.41888 -56578 1.57215 1.35435 -54231 2.72673 
75° 1.41453 -48458 1.57433 1.38417 -55061 2.70411 
1.41421 -44721 1.57607 1.39496 -55084 2.69164 
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TABLE II. Rotational eigenfunctions. 


Ve 








V,cm7! 


Oss 


28A 


2AS 


4ss 


4AA 


6SA 


6AS 


8AA 


8ss 


12AA 


16AA 





0.000 
377 
1.486 
3.283 
5.688 
8.601 
11.891 
15.425 
19.048 
22.614 
25.974 
28.991 
31.548 
33.543 
34.903 
35.579 
35.555 
34.840 
33.478 
31.535 
29.106 
26.303 
23.256 
20.109 
17.000 
14.079 
11.479 


2.330213 
2.303053 
2.222462 
2.091441 
1.915922 
1.705351 
1.472545 
1.232386 
0.999585 
-786517 
-601576 
-448620 
-327469 
-235073 
- 166842 
-117781 
-083247 
-059342 
.043014 
-031998 
-024675 
-019918 
-016950 
-015239 
-014406 
-014168 
-014300 
-014618 
-014967 
-015228 
-015324 


330308 
303147 
222551 
091524 
915995 
-705412 
472592 
-232418 
999601 
-786514 
-601554 
-448574 
-327392 
-234957 
- 166666 
2117523 
-082872 
-058798 
-042230 
-030873 
-023077 
-017682 
-013881 
-011118 
-009005 
-007276 
.005749 
-004306 
-002882 
-00 1447 
-000000 


One ee NN NN 


0.00000 
-00217 
.00431 
-00642 
-00850 
-01054 
-01259 
-01476 
01721 
-02024 
-02427 
.02983 
-03782 
-04936 
-06609 
-09024 
-12483 
-17376 
-24187 
-33473 
-45810 
-61702 
-81441 

1.04940 

1.31572 

1.60060 

1.88487 

2.14455 

2.35406 

2.49050 

2.53787 


—0.03222 
—.03171 
—.03013 
—.02749 
—.02387 
—.01943 
—.01426 
—.00854 
—.00235 

-00431 
-01160 
-02012 
-03050 
.04394 
.06213 
-08735 
-12273 
-17224 
.24076 
-33391 
45751 
-61659 
-81412 
1.04923 
1.31565 
1.60062 
1.88497 
2.14471 
2.35427 
2.49074 
2.53812 


0.00000 
.22343 
-43731 
-63155 
-79615 
-92225 

1.00383 

1.03933 

1.03255 
99219 
-93014 


—0.12917 
—.12548 
—.11357 
—.09364 
—.06654 


oo 
=e 
nm 
> CO 
oo 


-00755 
86630 
57647 
14540 
-60331 
-00000 


ee et ee DD ND et et et et ee 


0.00000 
-42040 
-82077 

1.18004 

1.48082 

1.70481 

1.84017 

1.88352 

1.84184 

1.72956 

1.56996 

1.38437 

1.19324 

1.01202 
-84933 
-71274 
.59960 
-50738 
-43204 
-36844 
-31156 
-25682 


0.0000 
-3677 
-7179 

1.0327 

1.2937 

1.4848 

1.5948 

1.6200 

1.5660 

1.4466 
1.2802 


—0.634 
—.610 
—.532 
—.403 
—.230 
—.024 

-202 
-428 
641 
-827 
.980 
1.098 
1.188 
1.256 
1.312 
1.362 
1.410 
1.556 
1.493 
1.510 
1.493 
1.422 
1.281 
1.057 
-748 
.369 
—.048 
—.456 
—.803 
— 1.036 
—1.118 


0.000 
.295 
-567 
.787 
.931 
977 
914 
-741 
473 
-137 

—.233 
—.605 
—.948 
— 1.243 
—1.470 
—1.626 
— 1.705 
— 1.700 
—1.592 
— 1.400 
— 1.097 
—.700 
.—.235 
.254 
.706 

1.054 

1.232 

1.199 
951 
525 
.000 


0.000 
481 
893 

1.175 

1.277 

1.170 
856 
.378 

—.190 
—.753 
—1.215 
— 1.504 
—1.579 
— 1.434 
— 1.097 
—.621 
—.068 
-487 
971 

1.314 

1.456 

1.363 

1.034 
513 
—.108 
—.700 
—1.127 
—1.276 
—1.100 
—.636 
.000 


20AA 


0.000 
-584 
1.042 
1.277 
1.220 
861 
-267 
—.428 
— 1.047 
— 1.420 
— 1.436 
— 1.079 
—.442 
.319 
1.010 
1.457 
1.549 
1.268 
-686 
—.049 





cm"! 


8.8865 


8.8866 


13.8952 


13.8962 


27.01 


27.54 


28.03 


28.81 


35.98 


45.4 


61.0 








r 


These solutions will be either symmetric or anti- 
symmetric at g=0 and at g=7/2, independ- 
ently; there are thus four symmetry types, 
which will be designated SS, SA, AS and AA. 
Several solutions of each type have been com- 
puted by the numerical method of Kimball and 
Shortley,? using an interval of 3°. In these 
solutions we have taken (1/27)(k/m)'/?=698 
em! and then have replaced >>W; by the 
simpler expression 


1878.8[ (1+0.154392y? —0.162276y4)!/? 


—0.996050!cm-'!, (43) 


which is not greatly different from it ; the process 
of solution uses only numerical values in any 
case, but the use of (43) avoids having to find 
the roots of (18) and (19) for 31 values of y. 
The calculated eigenfunctions and eigenvalues 
are presented in Table II and Fig. 3; the second 
column of Table II gives values of the potential 
energy term (43). The main quantum number 
characterizing each solution has been taken to 
be the number of nodes in a period 0 to 27; 
the symmetry type is shown as a subscript. 

A number of points related to these solutions 


‘Kimball and Shortley, Phys. Rev. 45, 815 (1934). 


are worth attention. The lowest level is a very 
narrow doublet, lying 8.89 cm above the 
potential curve at g=0, with a separation of 
about 1X10-* cm™, the Oss level being lower 
than 2s,; the eigenfunctions for these two states 
fall off very rapidly as ¢ increases from 0°; if 
g=45° is taken as the dividing line, the proba- 
bility of finding the molecule in the trans 
position is only 8X10~ for Ogg and 7X10~ 
for 23,4. The second level is again a narrow 
doublet, 7.54 cm™! above the potential curve at 
g=90°; 2as is 1010-4 cm— below 435; this is 
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Fic. 3. Rotational eigenfunctions and eigenvalues. 
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definitely a trans state, with a probability of 
only 510-4 for the molecule to be found in the 
cis form. 64g is cis and 6s, trans; these levels are 
still quite definitely vibrational in character, and 
have approximately three times the energy of 
the lowest cis and trans levels. It is to be expected 
that 4,4, and 8,4, will form doublets with these 
levels; since the latter two levels are of nearly 
the same energy, however, their would-be doub- 
lets, which are of the same symmetry type, can 
perturb each other. This is seen to occur, and 
neither 4,4, nor 8,4, can be characterized as 
either cis or trans, although they have far too 
little energy to execute classical rotations. All 
higher levels are actually rotating; 8ss barely 
exceeds the energy of the potential maximum, 
and hence the eigenfunction is largest in the 
perpendicular region near gy=7/4. The higher 
levels execute relatively free rotation. 

It is now possible to return to (40) and com- 
pute the first-order perturbation energy. It is 
easily found that for Ogg and 23,4, W'= —0.04 
cm, while for 24g and 4ss, W'= —0.01 cm. 
All the other rotational levels of the non- 
vibrating molecule will have values of this order, 
the maximum possible being —0.07 cm. Thus 
if the first-order perturbation energy is a reliable 
criterion, which of course may not be the case, 
the wave equation is for all practical purposes to 





d*Ye 
A 
ov 6. 


Ove. 2 
B—+ 
OV6 


h Vo 


instead of (31). This might be desirable if the 
calculation were to be carried to the second 
order, but for the present purposes it seemed 
better to use the obvious form (31) and get a 
first-order energy which one might hope would 
give some indication of the accuracy actually 
obtained. 


FREQUENCIES OF PLANE ForMs OF C; 


A treatment of the five-carbon chain which 
attempted as much as has been done for C, 
would be hopelessly difficult. Here there are 
seven vibrational and two rotational coordinates, 
and no partial separation into symmetric and 
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be considered separated in the form (29). For 
higher values of the vibrational quantum num- 
bers, particularly of mz and m;, the perturbation 
energy will be larger, but it can never exceed 
—0.07(n?-++-n+1) cm“ where x is the largest of 
the five n;. This is of the order of 10-? of the 
anharmonic term for most vibrations and is thus 
a negligible contribution to the absolute positions 
of the levels. The change in relative separations 
of the rotation levels is somewhat more sig- 
nificant, but scarcely serious for moderate values 
of the quantum numbers. 

The form of the pseudo-potential energy curve 
> W; will be different for each vibrational state. 
High values of 2; and mz tend to deepen the trans 
minimum and obliterate the cis one, while 1, 
and n; have the reverse effect; values of m3 of 8 
or more would give a shallow minimum near 60°. 
It is particularly important to note that every 
set of 1; consisting of four 0’s and one 1 gives 
minima at both the cis and the trans position, 
and hence probably has both cis and _ trans 
levels. There should accordingly be absorption 
bands and Raman lines corresponding to the 
five frequencies of each of the two plane struc- 
tures, that is, near 1106, 1100, 1044, 987, 974, 
879, 492, 385, 312 and 221 cm™. 

It is possible to make the first-order perturba- 
tion energy vanish, by taking 


n;+n,+1 A Op: 2 
eatttectt Agim) 
32 pi? OVs 





antisymmetric coordinates is possible. One would 
have to expand a ninth order determinant to find 
the seventh degree equation in which de- 
termined the frequencies as functions of the two 
angles and then solve numerically a two-dimen- 
sional wave equation to study the rotational 
motion. It seems likely from the work on Cs, 
however, that for C; also the plane structures 
will be of outstanding importance, and _ the 
vibration frequencies therefore have been calcu- 
lated for them. There are three such plane 
structures for C;, as may be seen in Fig. 4. 
Two of them, a W and a pentagon, have 4 
symmetry axis, and hence give normal modes of 
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vibration which are either symmetric or anti- 
symmetric; the third form, a dipper, has no 
helpful symmetry properties, and leads to an 
unfactored seventh degree determinant. Since 
the method of calculation is perfectly straight- 
forward, we give only the Hamiltonians and the 
expanded secular equations. The kinetic energy 
terms in the Hamiltonian will be useful for a 
possible subsequent calculation with a more 
general potential energy. 

The original coordinates are shown in Fig. 4. 
We make the substitutions 


6s+62=m1, 63— 62= no, 


54+61= 73, 64—61= 4, 
€3+€2=01, €3 — €2=02, 


Fic. 4. Coordinates for three plane forms of C;. s+ €1=03, €4— €1= 04. 


Then for the W-form 
T = (md? /60)(226;?+ 186163+2°!?¢19) —12-2"/?6 193 +1363? +27/2G391 — 2°/2o33+1991°? + 167173 
+11%3°) + (md?/636)(7062?-+ 266064 — 2°!2o94+91642—28-2"/2G 494 +41 942); (46) 
for the pentagon form 
T = (md2/540)(1986,? —306163-+300- 21/2691 —156- 21/2643 +8562 
— 80- 21/9637: 4-20-2"/?6¢392+395 91? — 4009193 + 131737) 


+ (3md2/1012)(70622+ 586264 — 52+ 2"26o%4+59642—36-2"26494+41 942); (47) 


and for the dipper 
T = (md? /22,860) (741461? + 246461¢2+47306163 — 3184616 4+7112-2"/26191 —3388-2"/6i73 
—1720-2"?6193+376662?+ 16006203+6586264+3556-2"/*G0n1 —1784-2"2Gons 
— 2660-22605 + 398565? —880c3¢4+ 2540-22639; —2200-2"/2¢ 393 —820-2"/2G m4 
+3991 64? —3556-2"/?6 491: 4+524- 2"? 6493+320-2"/?6494+11,30391? —40647173 
— 10,16091741+ 2453 94?+ 16007371+3935 9,7). 
In each case the potential energy is 
V = (d?/4) (koi? + hoo? +ko3?+hog?+2bn1?+bn3?+5n,”). 
The secular equations for the W-form are 
M—[(11/3)(k/m) + (29/3) (b/m) Js? +[ (29/9) (k®/m?*) + (247/9) (bk/m*) + (82/9) (b?/m?*) Jr? 
—[(41/3) (bk? /m*) + (70/3) (b?k/m*) JA+ 10b2k?/m*=0 





334 


and 


x3 —[ (13/3) (k/m) + (13/3) (b/m) Tr2-+[ (41/9) (k2/m2) + (161 /9) (bk /m2) JA— (53/3) (bk2/m3) =0; 


for the pentagon-form 
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(51) 


A4—[(11/3)(Rk/m) + (29/3) (b/m) JA8+[ (29/9) (k2/m?) +31bk/m? + (82/9) (b?/m?) Jr? 
—[(73/3) (bk? /m*) + (566/27) (b?k/m*) JA+10b2k?/m4=0 


and 


8—[ (13/3) (k/m) + (13/3) (b/m) ]A°+[ (41/9) (k2/m?) + (129/9)(bk/m?) Jr 


and for the dipper-form 


— (253/27) (bk?/m*) =0; 


7 —[8k/m+14b/m ]r\°+[(71/3) (k2/m?) + (928/9) (bk /m?) +51b?/m? ]rA> —[ (92/3) (k?/m?) 
+ (2486/9) (bk2/m) +344b°R /m®+ (1066/27) (b3/m?) JA4+ (1/81)[1189k4/m*+25,672bk3 /m! 
+67,303b7k?/m*+ 19,664b%k/m* ]r3 — (1/81)[10,550bk4/m>+ 68,156b2k?/m>+42,838b%k?/m> |r? 


+ (1/81)[24,079b?k4/m®+38,312b%k3/m® JA — (1270/9) (b%k*/m7) =0. 


This last equation comes from an unfactorable 
seven-rowed determinant which might appear 
impossibly laborious to expand; by suitable 
manipulation, however, it is possible to reduce 
the determinant one row at a time, always 
retaining somewhere a monomial term with 
which a row of zeros can easily be produced; 
the expansion in this way requires only three or 
four hours. 

When 6/k=0.1 and (1/27) (k/m)!/?=698 cm, 
the frequencies are as given in Table III. 


TABLE III. 








W-form 


1041 1092 1126 
1008 853 1062 
336 559 995 
213 143 908 
1115 1127 467 
946 981 343 
428 298 174 


Pentagon-form Dipper-form 











The first four frequencies for the W- and 
pentagon-forms belong to symmetrical, the last 
three to antisymmetrical, vibrations. 


COMPARISON WITH EXPERIMENT 


There are no data on the fine structure of 
absorption or of Raman lines which could be 
compared with the calculated torsional and low 


(54) 





rotational levels of the g-motion. Such a com- 
parison would in any case be difficult without 
knowledge of the selection rules, which we have 
not found, and without any theoretical study of 
the states with nonvanishing total angular mo- 
mentum. Hence it is only possible to seek 
necessarily rough correspondences between calcu- 
lated vibration frequencies and the observed 
spectra. The choice b/k=0.1 has been found to 
give the best results; Fig. 5 shows the calculated 
frequencies for C2, C3, Cy and C;, and such 
experimental data as are available for ethane, 
propane, butane and pentane. 

The C—C frequency for ethane is definitely 
located at 993 cm! by the Raman data of 
Bhagavantam® and. Lewis and Houston.’ The 
infrared line at 838 cm~ reported by Coblentz” 
cannot be understood, since it seems too low to 
be a C—H bending frequency. The Raman data 
on propane are from Bhagavantam and Kobhl- 
rausch and K6ppl," the infrared from Barto- 
lome;” an infrared line at 750 cm= which is 
tentatively identified as the harmonic of the 
Raman line at 373 cm~ has accordingly been 


8’ Bhagavantam, Ind. J. Phys. 6, 595 (1931). 

® Lewis and Houston, Phys. Rev. 44, 903 (1933). 

10 Coblentz, Carnegie Institute Publications, No. 35 
(1905). 

4 Kohlrausch and Képpl, Zeits. f. physik. Chemie B26, 
209 (1934). : 

'2 Bartolome, Zeits. f. physik. Chemie B23, 152 (1933). 
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Fic. 5. Calculated and observed frequencies. 


omitted from the figure. The Raman data on 
butane are due to Bhagavantam, Murray and 
Andrews and Kohlrausch and Képpl, while the 
infrared are from the old work of Coblentz. 
In the case of pentane, twelve Raman lines below 
1100 cm™! were reported by Ganeswan and 
Venkateswaran,“ a number of which have not 
been confirmed by the later work of Kohlrausch 
and K6éppl and are indicated as doubtful. No 
infrared data are available for pentane. The 
calculated deformation frequencies for the pen- 
tagon form are certainly worthless due to the 
repulsion between the terminal CH; groups. 

When allowance is made for the extreme 
approximations of the theory, and the un- 
certainty of the experimental data, there does 
seem to be a fairly satisfactory statistical agree- 
ment between the two. One may therefore feel 
that the frequencies calculated for still longer 
chains will be approximately correct and can be 
used to estimate specific heats and other thermo- 
dynamic properties without gross errors. It is 
intended to do this in the near future. 
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APPENDIX. THE KINETIC ENERGY OF A 
ROTATING MOLECULE 


Treatises on mechanics give the rotational 
energy of a rigid body in terms of moments of 
momentum about the principal axes. A more 


'’ Murray and Andrews, J. Chem. Phys. 1, 406 (1933). 
noes neswan and Venkateswaran, Ind. J. Phys. 4, 195 
). 
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general formula is derived here. Let the moments 
and products of inertia be 


A=m)>(y?+2?), B=m>o(x?+:2°), 
C=m>(x?+y?), D=mdyz, 
E=m)_ x2, F=m)> xy, 


»(55). 


and the instantaneous moments of momentum be 
M,=m)>(yz—zy), My=m>-(2x—42), 
M,=m>(éy—yx). 


It is supposed that the origin of the system is 
fixed at the center of gravity. Transform te 
rotating axes by the infinitesimal rotation 


x=X —e,Vi+e,Z1, 
y=e.Xt+ Y—.«,Z1, Y= —e,xt+y+e,2t, (57) 
2=—e,Xite,Vit+Z, Z=e,xt—e,yi+2, 


(56) 


X =x+e,yt—€,2t, 


where ¢ is at once the time and an infinitesimal 
parameter. Setting the moments of momentum 
in the new coordinate system equal to zero gives 
the conditions 


Ae,— Fe,—Ee,= —M,, 
— Fe,+Be,—De,= —M,, 
— Ee,—Deyt+Ce,= — M;. 


(58) 


It is then easily shown that 
Trot = (m/2) 2 (4? +9? +2") 
— (m/2)0(X?+ ¥+2?) 
= —(1/2)(e.M.zt+e,M,+e.M.) 
=(1/2)[(BC—D?)M/°+(AC-E’?)M/ 
+(AB—F*)M,?2+2(CF+DE)M.M, 
+2(BE+DF)M.M.+2(AD+CF)M,M.] 
xX [ABC—AD?—BE*?—CF?—2DEF}". (59) 
When only one product of inertia, say E, is not 
zero, the result reduces to 
Trot = (1/2) [M,?/B+(CM,?+2EM.M, 
+AM,’)/(AC—E?*)], (60) 


which is the form needed in this paper. 
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The Free Energy of Sulphur Dioxide 


A. R. Gorpon, Chemistry Department, University of Toronto 
(Received March 18, 1935) 


Values of (F°—E°%o)/T for sulphur dioxide are computed from spectroscopic data; the heat 
of formation of sulphur dioxide from diatomic sulphur and oxygen is calculated from the 
calorimetric heats of formation of sulphur dioxide and hydrogen sulphide, and Cross’ value for 
the energy of dissociation of hydrogen sulphide. Equilibrium constants for the reactions 


252(g) +O2.=SO2, 


are tabulated for the range 298.1 to 2800°K. 


SO + 302:=SO, 





F all the simpler gases the one which 

presented the greatest experimental diffi- 
culties in the determination of its free energy 
was probably sulphur dioxide. There was no 
appreciable dissociation at any attainable tem- 
perature, and the only values of its free energy 
of formation have been determined indirectly 
from other equilibria. * * The success that has 
attended the use of spectroscopic data in pre- 
dicting equilibrium constants for other reactions 
suggests that the same method might be used 
for sulphur dioxide, although in this case the 
data are not as complete as one could wish. 

For SO, the constants‘ are: A, B, C=12.3, 
73.2, 85.5X10-"; v(o) =1378—8v,, v(x) = 1150, 
6=520 cm. Ionescu gives no details of his 
rotational analysis, but his moments of inertia 
are consistent with Kohlrausch’s value® of the 
valence angle (119° 40’) and with Wierl’s value® 
for the S-O distance (1.37A). The resulting 
(F°—E°%o)/T are given in the second column of 
Table I. Since (F°— E°o)/T is known for sulphur 
vapor’ and for oxygen,® the only additional 
information that is needed for a calculation of 
equilibria involving SO, is the value of AE%y or 
of AHesg for the reaction 4S2(g)+O2.=SOs2. 


1 Lewis and Randall, J. Am. Chem. Soc. 40, 362 (1918). 
( . _ and Bichowsky, J. Am. Chem. Soc. 40, 368 
1918). 

3’ Ferguson, J. Am. Chem. Soc. 40, 1626, 1900 (1918). 

4Tonescu, Comptes rendus 196, 1476 (1933); Smyth, 
Phys. Rev. 44, 690 (1933). The universal constants used 
in the calculations of this paper are: h=6.554X10-?’, 
k=1.372X10-®, c=2.9986X10', R=1.9869, additive 
constant for the translational free energy = — 7.267. 

5 Kohlrausch, Physik. Zeits. 33, 165 (1932). 

6 Wierl, Ann. d. Physik 8, 521 (1931). 

7 Montgomery and Kassel, J. Chem. Phys. 2, 417 (1934). 
=_ and Walker, J. Am. Chem. Soc. 55, 172 
1933). 


TaBLe I. Ki=(Pso:)/(Ps2)'!*(Po,) i 
K2=(Pgo:)/(Pso)(Po2)'. 








—(F° —E%)/ —(F°—E%)/ 
T ‘ T 


log K2| T°K 


$4.53 900 
54.16 | 1000 
45.84 | 1200 
39.60 | 1400 
30.85 | 1600 
25.01 | 2000 
20.84 | 2400 
17.72 | 2800 


T°K 


298.1 
300 50.62 
350 $1.93 
400 53.10 
500 55.11 
600 $6.83 
700 58.33 
800 59.67 


log Ki log Kz 





50.56 -59.65 
59.25 
50.24 
43.48 
34.02 
27.70 
23.19 
19.80 


60.89 
62.01 13.34 
64.00 10.42 
65.73 8.33 
67.27 6.77 
69.92 4.59 
72.15 3.13 
74.06 2.09 


15.28 











Eckman and Rossini? have measured the heat of 
combustion of rhombic sulphur (AJ/293 = — 70,940 
+50 cal.) but thermal data for the reaction 
S(rh) =4S2(g) are unknown. These can be ob- 
tained, however, with sufficient accuracy for the 
purpose in hand from the recent calculation of 
Cross! for the hydrogen sulphide dissociation 
based on his rotational analysis of the hydrogen 
sulphide spectrum,!! and the determination by 
Zeumer and Roth” of the heat of formation of 
hydrogen sulphide from rhombic sulphur and 
hydrogen (A//29;= —4800+150 cal.). Cross finds 
for the reaction H2+3S2(g)=H2S, AE, = —19,- 
620 cal.; from this, and the spectroscopic free 
energies of hydrogen sulphide,!® sulphur!® and 
hydrogen,!4 AF%o93.1=—17,570 cal. and_ since 
S°o98.1 for HeS, Ss and He is 49.15, 54.42 and 
31.23, respectively, Als is — 20,340 cal. Com- 


®Eckman and Rossini, Bur. Standards J. Research 3, 
597 (1929). 

10 Cross, J. Chem. Phys. 3, 168 (1935). 

1 Cross, Phys. Rev. 46, 536 (1934); 47, 7 (1935). 

12Zeumer and Roth, Zeits. f. Elektrochemie 40, /// 
(1934). 

13—n the calculations of this paper, Montgomery ana 
Kassel’s values of —(F°—E%)/T (reference 7) have been 
increased by 0.55 cal./degree as recommended by Cross 
(reference 10). 

14 Giauque, J. Am. Chem. Soc. 52, 4816 (1930). 
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bining this value with Zeumer and Roth’s result, 
it follows that A/Jo9, for the reaction S(rh) 
=3S2(g) is 15,540 cal., and consequently AJJ29s 
for the reaction 3S2(g) -+O2=SOs, is — 86,480 cal. 
Since S°9s.1 for oxygen® and for sulphur dioxide 
is 49.02 and 59.02, respectively, AF°29s.1 for the 
formation of SO, from diatomic sulphur and 
oxygen is — 81,350 cal. and AE% is — 85,860 cal. 

There is an alternative method of obtaining 
the heat of vaporization of diatomic from 
rhombic sulphur similar to that used by Lewis 
and Randall.!® For the four reactions 

H.S=H2+3S,(g), AF°o98.1= 17,570 cal., 
2HI+S(rh) = HS+I2(s), ‘= — 8470 cal.,'® 

H.+1.(g) =2HI, ‘* = — 3690 cal.,!? 

I,(s) = I2(g), «= 4620 cal.!8 
Hence for the reaction S(rh)=3S.(g), AF°o98.1 
= 10,030 cal., and, since the third law entropy 
of rhombic sulphur is 7.69 cal./degree,'!® AT/o9s 
for this reaction is 15,850 cal., which agrees with 
the previous result (15,540 cal.) probably as 
closely as can be expected. 

In the third column of Table I are given 
the values? of log Ki=log (Pso,)/(Ps,)*(Po,), 
computed from the spectroscopic values of 
(F°—E%)/T for the constituents, on the as- 
sumption that AE%, for the reaction is — 85,860 
cal. In the last column are the equilibrium 
constants for the “sulphur monoxide’’ dissocia- 
tion of SO2, log K2=log (Pgo,)/(Pso)(Po,)?; these 
were obtained by means of Montgomery and 
Kassel’s table of (F°—E%)/T for SO, on the 
assumption that AE% for the reaction SO+40, 
=SO. is —79,300 cal. It is evident that the 
equilibrium constant for the reaction 45S,(g) 
+SO,=2S0O, is given by the relation log (Pgo)?/ 
(Ps,)'(Pso,) =log Ki—2 log Ks. 

Of the equilibria involving SO, which have 
been studied experimentally, Randall and Bi- 
chowsky’s measurements? for the reaction SO. 
+3H.=2H20+H2S and Ferguson’s* for the re- 


Ne Lewis and Randall, Thermodynamics, 1923, p. 540 
et seq. 

'® Pollitzer, Zeits. f. anorg. Chemie 64, 121 (1909). 

" Zeise, Zeits. f. Elektrochemie 40, 885 (1934). 

's Baxter, Hickey and Holmes, J. Am. Chem. Soc. 29, 
127 (1907) ; Giauque, J. Am. Chem. Soc. 53, 507 (1931). 

Int. Crit. Tab. 5, 89. 

**A preliminary table of log K; is given in a recent 
review article by Zeise, Zeits. f. Elektrochemie 40, 885 
(1934). The values given there differ slightly from those of 
this paper, since at that time Cross’ calculations for 
hydrogen sulphide had not appeared. 
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action SO.+2CO=3S,(g)+2CO. can be com- 
pared with the results of calculation. For the 
former, from Rossini’s value®! for the heat of 
formation of steam, AE%, for the reaction is 
— 47,980 cal.; from this and the known free 
energies of hydrogen,'* steam” and hydrogen 
sulphide,"° log (Pu,0)?(Pu,s)/(Pso,)(Pu,)* for 
1160°, 1362°, 1473° and 1645°K is 5.58, 4.06, 
3.40 and 2.54, respectively. Randall and Bi- 
chowsky found for the same four temperatures 
5.93, 4.32, 3.54 and 2.56. Calculation shows that 
the amount of sulphur monoxide present under 
the conditions of their experiments is quite 
negligible, and also (from the data of Mont- 
gomery and Kassel) that the dissociation of 
diatomic into monatomic sulphur can likewise 
be ignored. There is thus a discrepancy of about 
1.6 cal./degree in Rln K at the lowest of their 
temperatures, this difference decreasing rapidly 
as the temperature rises. This improvement in 
the agreement is perhaps significant, since one 
would expect caeleris paribus that equilibrium 
would be the more readily established the higher 
the temperature. It should be noted in passing 
that the revised free energy equation of East- 
man?’ for SOs combined with the spectroscopic 
free energies of formation of HzO and H.S gives 
for log K at the same four temperatures, 6.50, 
4.97, 4.30 and 3.43—definitely in worse agree- 
ment with Randall and Bichowsky’s results than 
are those obtained above. 

For the second reaction, studied by Ferguson, 
from Rossini’s heat of combustion of carbon 
monoxide*‘ and the spectroscopic free energies,”* 
log . (Pso,)(Pco)?/(Ps,)*(Pco,)? for 1275° and 
1460°K is —3.07 and —2.02, respectively; 
Ferguson found —3.56 and —2.38. The dis- 
agreement between the results of calculation and 
experiment is thus definitely greater than the 
apparent uncertainty in the experimental data 
(+0.2 in log K) on the one hand and in the 
calculated thermodynamic quantities on the 
other, although here again the discrepancy 
decreases with rise in temperature. 


2) Rossini, Bur. Standards J. Research 6, 1 (1931). 

22 Gordon, J. Chem. Phys. 2, 65, 549 (1934). 

23 Eastman, Bur. Mines Infor. Circ. No. 6454 (1931). 

24 Rossini, Bur. Standards J. Research 6, 37 (1931). 

25 For CO, Clayton and Giauque, J. Am. Chem. Soc. 55, - 
iChh net for COs, Kassel, J. Am. Chem. Soc. 56, 1838 
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The Mercury Photosensitized Polymerization of Acetylene and Acetylene-d, 


JosePH C. JUNGERS! AND HuGu S. Taytor, Frick Chemical Laboratory of Princeton University 
(Received March 23, 1935) 


The mercury photosensitized polymerization of acetylene 
and acetylene-d, has been studied. The rate of reaction is 
30 percent greater with acetylene than with the deutero- 
compound, over the whole pressure range investigated. 
At low pressures, the rate varies as the acetylene pressures 
but in the pressure range 0.7-4 cm the initial rate of 
polymerization is practically independent of the acetylene 


pressure and determined by the intensity of absorbed 
resonance radiation. The quantum yield is approximately 
6.5 for acetylene and 5 for acetylene-d, lower than the 
yield for acetylene at the same temperature in the photo- 
chemical reaction (A~2000A). The results are also com- 
pared with corresponding studies of the polymerization 
process under the influence of a-particles. 





E have applied the technique developed 

by us in our studies of the mercury 
photosensitized decomposition of the deutero- 
ammonias’ to a comparison of the polymerization 
of acetylene and acetylene-d2 under the influence 
of excited mercury. Bates and Taylor® indicated 
qualitatively that acetylene could be polymerized 
to solid cuprene under such conditions. The 
straight photochemical polymerization of acety- 
lene has already been studied by Lind and 
Livingston who have determined the influence 
of pressure and temperature and the quantum 
yield, this latter being compared with the ion 
yield in the parallel determinations of polymeri- 
zation under the influence of a-particles.> One 
of us (J.C.J.) has already, in collaboration with 
Lind, compared the rates of polymerization of 
acetylene and acetylene-d. under the influence 
of a-particles and found that these rates were 
identical within the error of experiment.® It was 
of interest to ascertain whether this identity of 
reaction rates extended to the reactions induced 
by light. 


APPARATUS 


The experiments were carried out with the 
apparatus already described by us in our record 
of the work on the deuteroammonias.? The 
quartz reaction-vessel there employed was very 
suitable also in the present work owing to the 


1C.R.B. Fellow from Louvain University, 1934-35. 
2 Jungers and Taylor, J. Chem. Phys. 2, 373 (1934). 
3 Bates and Taylor, J. Am. Chem. Soc. ‘49, 2444 (1927). 
4 Lind and Livingston, J. Am. Chem. Soc. 54, 94 (1932). 
5 Mund and Koch, Bull. Soc. Chim. Belg. 34, 125, 241 
: eee Lind and Bardwell, J. Am. Chem. Soc. 48, 1556 
(1926). 

6 Lind, Jungers and Shiflett, Phys. Rev. 46, 825 (1934). 


device employed for increasing the sensitivity of 
the pressure measurements which were made 
after compression of the reacting gases to a small 
volume. This permitted us to study accurately 
small amounts of polymerization at relatively 
low pressures and thus to avoid the difficulties 
attendant on photochemical experiments in 
which opaque deposits of product form on the 
reaction vessel. As source of resonance radiation 
we again employed the Hanovia quartz mercury 
vapor, rare gas discharge tube, Sc—2537, oper- 
ating at 5000 volts from the 110-volt a.c. lighting 
circuit through a 30-milliampere transformer. 
The advantages of this type of tube have already 
been described.? We used the radiaticn trans- 
mitted through a quartz cell containing 25 
percent acetic acid, 1.5 cm in thickness which 
eliminated the short wavelength light that is 
photochemically active, but which transmits the 
resonance radiation. Between each experiment, 
the quartz reaction vessel was cleaned by heating 
it in the presence of oxygen. 


MATERIALS 


The samples of acetylene and acetylene-d2 were 
prepared by interaction of the vapors of water 
and of pure deuterium oxide (d2;2°= 1.1079) with 
samples of pure calcium carbide. This carbide 
was placed at our disposal by the Union Carbide 
and Carbon Research Laboratories, Inc., through 
the courtesy of Dr. C. E. MacQuigg. The sample 
which the donor believed to be quite pure was 
entirely different in outward appearance from 
commercial products. It was nut-brown in color 
and was preserved in a moisture-proof can. For 
the preparations, the inner cores of freshly 
broken lumps were chosen, introduced into the 
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preparation vessel, heated and exhausted in 
vacuum at 400°C for 24 hours before reaction 
with the water vapors. The products were con- 
densed in liquid air and the middle fractions only 
used for experimental work. We cite the data of 
Glockler and Davis,’ on the Raman spectrum 
of our acetylene-d; preparation as evidence of 
its distinctive characteristics as compared with 
acetylene and acetylene-d). 


EXPERIMENTAL RESULTS 


The rates of polymerization of the two 
acetylenes, at pressures of 0.1, 0.25, 0.7 and 4 
cm, are shown graphically in Fig. 1. The meas- 
urements were made at room temperatures 
(ca. 25°C) with the gases confined over mercury 
and hence saturated with mercury vapor at that 
temperature. It is immediately evident that, 
over the pressure range employed, the rate of 
acetylene polymerization is always greater than 
that of the acetylene-d.. Computation shows 
that, at each pressure, the rate, in the case of 
the light acetylene, is 30 percent greater than 
that of the heavy acetylene. The curves at the 
lowest pressure represent approximately 90 per- 
cent total reaction. At the higher pressures the 
change is not followed so far. It is also evident 
that, at the lower pressures, the initial rate of 
polymerization is strongly dependent on the 
pressure, but the data at the higher pressures 
show that the rate has become practically 
independent of the acetylene pressure. This 
pressure dependence is undoubtedly governed 
by the effectiveness of the acetylenes in taking 
up the energy of the excited mercury and 
reveals that there is no essential difference 
between the two acetylenes in this respect. In 
the higher pressure range the rates are therefore 
determined by the concentration of excited 
mercury 

—d(C2He) /dt= kL Hg’ ], 
—d(C,De) /dt = kL Hg’ ]. 


The available intensity of resonance radiation 
was ascertained by measurement of the rate of 
reaction of an electrolytic gas mixture (2H2+Oz) 


under identical operating conditions. It was 
found that the rate of polymerization of acetylene 


‘Glockler and Davis, Phys. Rev. 46, 535 (1934). 
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Fic. 1. Velocity of polymerization of acetylene and 
acetylene-ds at (A) 0.11, (B) 0.25, (C) 0.7 and (D) 4 cm 
initial pressure. 


in the high pressure range was approximately 
30 percent greater than that of the hydrogen- 
oxygen reaction under the same conditions of 
illumination and mercury pressure. Since the 
quantum yield in this latter reaction is found to 
be less than 3H,O2 per quantum® or approxi- 
mately 5 molecules of hydrogen plus oxygen 
disappearing, we assign a quantum yield of not 
more than 6.5 to the photosensitized polymeri- 
zation of acetylene and a 30 percent lower yield, 
or approximately 5 molecules per quantum for 
heavy acetylene. 


DISCUSSION 


The experiments which we have performed 
serve still further to differentiate the polymeri- 
zation of the acetylenes under the action of ions 


8 Marshall, J. Phys. Chem. 30, 34, 1078, 1634 (1926); 
J. Am. Chem: Soc. 45, 4460 (1932). Frankenburger and 
Klinkhardt, Zeits. f. physik. Chemie 8, 138 (1930); 15, 
421 (1932). 
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and of light. Whereas with the former (a- 
particles), there is no difference in the rates of 
polymerization with the heavy and light mole- 
cules, in the present work the heavy acetylene 
polymerizes at a markedly slower rate under the 
same intensity of excited mercury. There is a 
progression, also, in the yield of the respective 
processes. Under the influence of {-particles,° 
26 molecules of acetylene polymerize per ion 
pair. With a-particles the rate of disappearance 
of the two acetylenes is 18-20 molecules per ion 
pair’? practically independent of temperature. 
Lind and Livingston found a value for the 
photochemical polymerization (A~2000A) of 9.5 
molecules per quantum at 25°C and values 
varying from 7 to 11 per quantum in the temper- 
ature range of 7-40°C, with the temperature 
coefficient of 1.25 per 10° increase in temperature 
responsible for this variation, the absorption 
coefficient being temperature independent within 
the limits of accuracy of the measurements. 
With the smaller quantum of light energy 
involved in the photosensitized polymerization 
(A= 2537A), the quantum yield falls still further 
to 6.5 molecules per quantum at 25°C. 

All these data point to the conclusion that the 
polymerization process in the case of the acety- 
lenes rapidly dissipates the energy which is 
employed to initiate the process. With agents 
such as a-particles, B-particles and cathode rays 
the chain mechanism is most pronounced and 
the velocities of polymerization identical in the 
case of light and heavy acetylenes. In the case 
of acetylene radiated by a-particles it is known 
that the only observable effect of temperature is 







































— and Jungers, Bull. Soc. Chim. Belg. 40, 158 
(1931). 

10 An explanation for this lower yield with a-particles 
has already been put forward by Mund and Jungers, 
reference 9. Recombination of ions, due to their greater 
concentration in the path of the a-particle, may account 
for the lower value as compared with @-particles. 
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slightly negative (2 percent in 100°) and this 
has been explained on the basis of decreasing 
capture of molecules by the charged particle 
with increase in temperature." We conclude that 
the electrostatic attraction of the charged par- 
ticles produced by the radiation determines the 
magnitude of the ion yield and serves in place 
of the activation energy of collision necessary 
where such attractive forces are absent. With 
less energetic agents such as light, temperature 
activation is a markedly contributory factor, the 
evidence indicating that the heavy acetylene 
requires greater thermal activation for the chain 
process. In this respect, there is a possibility 
that the exchange of energy between excited 
mercury and acetylene-d2 is not so efficient as 
with acetylene. The absence of pressure effect 
at the higher pressures makes this, however, 
appear improbable. The effect of temperature 
as an auxiliary agent cannot be followed into 
the higher temperature ranges since it is known 
that in the region of 250°C the product is no 
longer cuprene, a complex of products, including 
benzene, being formed.” In the thermal poly- 
merization process, which requires an activation 
energy of the order of 40,000 calories,'* cuprene 
is also not produced but a complex mixture of 
products. All these polymerization processes as 
well as those of methyl acetylene and allene 
appear to involve only short chains, in contrast 
to the polymerization of styrene and vinyl 
acetate in which the chain length may reach 
thousands of molecules per quantum absorbed.” 


1 Mund, Ann. Soc. Sci. Brussels B54, 30 (1934). 

12 Kato, Bull. Inst. Phys. Chem. Res. (Tokyo) 10, 343 
(1931); Kemula and Mrazek, Zeits. f. physik. Chemie 
B23, 358 (1933); Livingston and Shiflett, J. Phys. Chem. 
38, 377 (1934). 

13 Pease, J. Am. Chem. Soc. 51, 3470 (1929). 

14 Lind and Livingston, J. Am. Chem. Soc. 55, 1036 
(1933). 

16 Taylor and Vernon, J. Am. Chem. Soc. 53, 2527 
(1931). 
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The Ethane-Ethylene-Hydrogen Equilibrium* 


HiLton A. SMITH AND WILLIAM E, VAUGHAN, Chemical Laboratory of Harvard University 
(Received March 18, 1935) 


An attempt has been made to correlate several data 
on this system by a statistical calculation of the equi- 
librium constants in the temperature range of 400 to 
700°C. The data are in good agreement but the calculated 
values of the K’s are smaller than the experimental by a 


factor of two. It is thought that the principal reason for 
this discrepancy lies in the calculation of the entropy of 
free rotation in ethane, although other factors may 
contribute. 





HERE are in existence two fairly recent 

determinations of the equilibrium constants 
of the ethane-ethylene-hydrogen system. The 
first is that of Pease and Durgan,' who, using a 
static system and approaching the equilibrium 
from both sides, found K’s at 600°, 650° and 
700°C. The other, which is the work of Frey and 
Huppke,’ was carried out in a flow system by 
catalytic dehydrogenation of ethane at one 
atmosphere pressure in the temperature range 
from 400-500°C. The results of these two re- 
searches are in good agreement.’ The purpose of 
this paper is to attempt to correlate by means of 
statistical calculations the figures obtained from 
these experiments with the recent accurate 
measurement of the heat of hydrogenation of 
ethylene at 82°C as determined by Kistiakowsky 
et al. 

For the present calculation the values of the 
entropy (excluding nuclear spin) and of the heat 
content of hydrogen at the desired temperatures 
have been taken by direct graphical interpolation 
of the data of Davis and Johnston.‘ 

There are in the literature several recent 
discussions of the computation of the entropy of 
polyatomic molecules such as ethylene and 
ethane.® It suffices for the present paper to give 
briefly the specific values and assumptions 
actually used. 


* Presented at the N. Y. meeting of the Am. Chem. 
Soc., 1935. 

‘ Pease and Durgan, J. Am. Chem. Soc. 50, 2715 (1928). 

* Frey and Huppke, Ind. Eng. Chem. 25, 54 (1933). 

* Kistiakowsky, Romeyn, Ruhoff, Smith and Vaughan, 
J. Am. Chem. Soc. 57, 65 (1935). 
ass and Johnston, J. Am. Chem. Soc. 56, 1045 

°See for example (a) Kassel, J. Am. Chem. Soc. 55, 
1361 (1933); (b) Mayer, Brunauer and Mayer, J. Am. 
Chem. Soc. 55, 37 (1933); (c) Halford, J. Chem. Phys. 2, 
—— This latter paper gives the formulae quite 
explicitly, 


The calculation of the translational entropies 
of ethane and ethylene requires no discussion; 
the values obtained are given in Table II. The 
estimation of the magnitude of the vibrational 
entropies is, however, a somewhat more difficult 
matter. Neither of these molecules has been 
completely analyzed. In the case of ethylene 
Kassel** has lumped Mecke’s® eleven frequencies 
in the following manner: 950 cm (triple) ; 
1100 (single) ; 1400 (double) ; 1620 (single); and 
3090 (quadruple). Kassel assumes for the fre- 
quency of the torsion about the double bond the 
value of 950 cm@'. Penney,’ however, gives 
750 cm (a figure quoted to him verbally by 
Sutherland). In view of the uncertainty of this 
frequency it was decided to select 790 cm™ 
which gives the best agreement with the accurate 
experimental values of Cp for ethylene as de- 
termined by Eucken and Parts,* as shown by 
Table Ia. It is seen that the figures deviate but 
slightly and in a random manner over quite a 


TABLE I. 








(a) 
Ethylene: vibrational heat capacity (12 degrees of freedom) 
T 143 160 196 245 280 373.5 464 
TCyid.(cale.) 0.092 0.183 0.540 1.347 2.078 4.164 5.914 


Cyib(E.&P.)* 0.08 0.21 0.53 1.30 2.03 3.96 6.20 
Deviation +0.01 --0.03 +001 +0.05 +0.05 +0.20 —0.29 








* 8R/2 has been deducted from the Cp values of Eucken and Parts. 


(b) 
Ethane: vibrational heat capacity (17 degrees of freedom) 
143 160 196 245 280 


373.5 





=C vib. (cale.) 0.242 0.434 1.043 2.224 3.211 6.068 
Cyib(E.&P.)* 0.32 0.50 1.16 2.13 3.18 5.89 
Deviation —0.08 —0.07 —0.12 +0.09 +.03 +0.18 





*9R/2 has been deducted from the Cp values of Eucken and Parts. 








® Mecke, Zeits. f. physik. Chemie B17, 1 (1932). 

7 Penney, Proc. Roy. Soc. A144, 166 (1934). 

8 Eucken and Parts, Zeits. f. physik. Chemie B20, 190 
(1933). 
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TABLE II. 
i 298° 355 673 773 873 973 
(a) Entropy of ethylene 
Translational 35.968 36.833 40.018 40.706 41.307 41.848 
Vibrational 0.671 1.233 5.761 7.305 8.831 10.320 
Rotational 15.909 16.439 18.349 18.759 19.119 19.439 
Total 52.548 §4.505 64.128 66.770 69.257 71.607 
(b) Entropy of ethane 4 
Translational 36.175 37.040 40.225 40.913 41.514 42.055 
Vibrational 1.149 1.950 8.106 10.220 12.320 14.390 
Rotational 18.903 19.602 22.146 22.694 23.180 23.611 
Total 56.227 58.592 70.477 73.827 77.014 80.056 
(c) Entropy of hydrogen 
Total 31.230 32.420 36.905 37.873 38.473 39.517 
(d) 
AS for H2+CoH«—-CoHe —27.551 — 28.333 — 30.556 — 30.816 — 30.986 —31.068 
(e) Heat content 
He 2023 2417 4642 §348 6057 6775 
CoHa 2534 3171 8062 9976 12,025 14,197 
CoHe 2939 3713 9768 12,195 14,822 17,625 
(f) AH for CoHs+H2—-C2He —32,567 --32,824+50 —33,885 —34,078 —34,209 — 34,296 
(g) AF for CoHs+H2—-CeHe — 24,358 —22,767 — 13,322 — 10,255 —7157 — 4069 
(h) K calc. 7.08( 1019) 1.01(10'4) 2.10( 104) 7.89( 102) 61.7 8.19 
(i) 1/K (calc.)* 0.000048 0.0013 0.016 0.122 
(j) K of CoHs—-CeHs+He (obs.) 0.000152 0.00322 0.0311 0.201 
0.000083 0.0024 0.031 0.17 


(k) K of CoHe—-C2Hi+He (recalc.)t 








* or K for the reaction written CoHse—CoHs+H2. 
tas done by Kistiakowsky et al.; see reference 3. 


large temperature range. Ethane has been treated 
in similar fashion using the assignments of 
Eucken and Parts; namely, 712 cm (double); 
826 (double); 990 (single); 1460 double; 1465 
(double); 1499 (double); and 2975 (sixfold); a 
total of seventeen vibrational degrees of freedom. 
The results obtained with these assignments are 
given in Table Ib. Use of these two sets of fre- 
quencies gives the values for S,i,, to be found in 
Table II. The entropy of nuclear spin has been 
omitted in all calculations. 

After a consideration of the various inter- 
atomic distances given by different authors, the 
following values were chosen for the calculation 
of the moments of inertia of the two hydro- 
carbons: C—C, 1.53A; C=C, 1.36; C—H, 1.08. 
The tetrahedral angle was assumed in the case of 
ethane, and the angle between the two carbon- 
hydrogen linkages on either end of the ethylene 
molecule was taken to be 110 degrees. This latter 
value is probably incorrect, but is sufficiently 
accurate for the present work. The moments of 
inertia found for ethylene were 5.205 (10-*°), 
29.53 (10-*°) and 34.71 (10-*°). Those used in 
the case of ethane were 10.27 (10-*°) and 
35.98 (10-*°) (double). A symmetry number of 


four was used in obtaining the values for Syo¢. of 
ethylene as given in Table IIa. 

The free rotation of the two methyl groups in 
ethane complicates the statistical computation of 
its entropy of rotation. Eucken and Weigert’ 
found their values of Cp for ethane to be in 
accord with the quantum-mechanical treatment 
of this problem as analyzed by Teller and 
Weigert.!° The latter showed that this contribu- 
tion to the rotational heat capacity becomes 
classical at 298°K, going through a slight 
“hump” in approaching this temperature. If one 
calculates the rotational entropy assuming free 
rotation to lowest temperatures (using two pairs 
of moments of inertia; namely, 5.14 (10-*°) and 
35.98 (10-°), and a symmetry number of 18), 
he obtains the data given in Table IIb. 

That the assumption of free rotation to lowest 
temperatures is permissible may be shown in 
the following manner. If C (of internal rotation) 
as obtained from the paper of Teller and Weigert 
is plotted against In 7, graphical integration 


®Eucken and Weigert, Zeits. f. physik. Chemie B24, 
277 (1934). . 

10 Teller and Weigert, Nachr. d. Ges. d. Wiss. zu 
Géttingen, 218 (1933), 
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gives the contribution to the rotational entropy 
at any temperature. This area was found to 
correspond to 7.241 e. u. at 298°K. The classical 
calculation of this portion of S,or. (setting = 9) 
gives the value of 7.160 e. u. These values are 
in reasonable agreement. 

The heat contents of Table IIc were arrived 
at by use of the vibrational frequencies already 
mentioned, together with the assumption that 
rotation has already become classical at these 
temperatures. The value of A// for the reaction 


H.+ C.H,—-C2H, 


on which the final data are based is —32,824 
+50 cal. at 355°K.* The summary of the calcula- 
tions is shown in the last few lines of the listing. 
It is immediately apparent that the three 
researches are essentially in agreement. How- 
ever, the calculated values of 1/K are fairly 
uniformly one-half as large as the experimental 
figures for the equilibrium constant. This dis- 
crepancy presents a real problem. To decrease the 
Kate. (line h) by a factor of two would necessitate 
a change of A/T beyond probability. At 298°K 
such an alteration would involve a change in AS 
of 1.376 e. u., and so would decrease the original 
difference in the entropy of ethane and ethylene 
from 3.7 to 2.3 e. u. This latter figure agrees 
well with the empirical rule of Parks and Huff- 
man" that the entropy difference between a 
saturated hydrocarbon and its corresponding 
olefine is 2.7 e. u. at 298°K. 
An attempt to localize this difference of 1.38 
e. u. in some part of the entropy calculation 
results in the elimination of Strans.. It likewise 
seems improbable that there is any large error in 
the estimation of S,;,. for the two hydrocarbons; 
indeed the contributions from this type of motion 
are small at 298°K and are undoubtedly of 
correct relative magnitudes. It is possible, how- 
ever, that at higher temperatures this calculation 


_" Parks and Huffman, Free Energies of Some Organic 
Compounds, Chem, Cat. Co., N. Y., 1932, p. 79. 
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may be in error. However, the inconsistency may 
more probably be attributed to S,o¢.. The con- 
figuration given ethylene in these calculations is 
certainly open to question. In fact Dieke 
and Kistiakowsky” from an analysis of band 
spectra of formaldehyde found 2.941 (10-*°) 
as the moment of inertia of the top. The value 
2 (2.941)10-* is larger than the one here calcu- 
lated, (5.205)10-*, as might be expected because 
of the attraction of the oxygen for the two 
hydrogen atoms. If one still uses 1.08A for the 
C—H distance, this experimental moment corre- 
sponds to an angle of 121° between the C—H 
bonds. But in the case of ethylene this spreading 
effect should be less, and hence the value 
5.88 (10-*°) might be taken as an upper limit. 
In fact Badger" has determined the moments of 
inertia of ethylene from spectroscopic data, and 
finds for the three moments the values of 33.2, 
27.5 and 5.7 (10-*). These give a product of 
5200 (10-'*°), while the value actually used was 
5330 (10-°). This discrepancy is quite neg- 
ligible. 

The model of ethane is, perhaps, less doubtful ; 
the interatomic distances are quite well estab- 
lished, and the tetrahedral angle must be a good 
approximation. The former factor should be 
sufficiently accurate to allow no errors greater 
than a few tenths of an entropy unit, and re- 
garding the latter, elimination of the discrepancy 
would require the assignment of fantastic values. 
Were the entire fault in the moments of inertia, 
the product of the four principal moments 
(assuming classical free rotation) would have to 
be reduced from 34,210 (107°) to 24,660 (10-') ; 
this is obviously too radical. 

Thus, although other factors may be con- 
tributory, most of the uncertainty would appear 
to lie in the treatment of the problem of the 
free rotation of the two methyl groups of 
ethane. 


12 Dieke and Kistiakowsky, Phys. Rev. 45, 4 (1934). 
18 Badger, Phys. Rev. 45, 648 (1934). 
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The thermal decomposition of n-propyl nitrite is a 
homogeneous first order reaction in the temperature range 
170 to 210°C. The mechanism of the decomposition is 
analogous to those of the other alkyl nitrites, viz., 


C;H;ONO—NO+3C:H;CHO+3C;H;OH. 
The rate of reaction is given by 


K =2.75 X10"e—37650/RT sec—1, 
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The rates and activation energies of the various nitrite 
decompositions are compared. It is concluded that the 
increase in rate with the higher members of the homologous 
series is probably due to the participation of extra vibra- 
tional degrees of freedom. An estimate of the strength of 
the —O—N-— bond is made. 








INTRODUCTION 


HE investigation of the unimolecular de- 

composition of a series of chemically related 
compounds is of considerable interest. It has 
recently been shown in the first two papers of 
this series! that the alkyl nitrites constitute an 
ideal series of compounds from this point of 
view. Methyl and ethyl nitrites were shown to 
decompose homogeneously and unimolecularly. 
The decompositions are simple and appear to 
involve the simple splitting of the —O—NO 
bond. 

Since submitting this paper, our previous 
conclusions regarding the mechanism of the 
reaction have been completely confirmed by a 
paper by Rice and Radowskas.? By the free 
radical method they have shown that, as we 
suggested, the first step in the decomposition of 
ethyl nitrite is 


C,H;NO:=C;H;0+ NO. 


Furthermore, their value of the activation energy 
for the split into radicals agrees with ours for 
the overall reaction, and hence confirms our 
conclusion that this split represents the rate 
determining step. 

In consequence it is possible to make com- 
parisons between the activation energies and 
reaction rates of the different members of the 
series. It is therefore desirable to investigate as 


1Steacie and Shaw, (I) Proc. Roy. Soc. A146, 388 
(1934); (II) J. Chem. Phys. 2, 345 (1934). 
2 Rice and Radowskas, J. Am. Chem. Soc. 56, 214 (1935). 
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many of the alkyl nitrites as possible. The 
present paper deals with the third straight chain 
member of the series, viz., n-propyl nitrite. 


EXPERIMENTAL 


Reaction velocities were measured as before 
by observing the rate of pressure change in a 
system at constant volume. The reaction vessels 
were Pyrex bulbs having a capacity of about 
125 cc. The apparatus was similar to that used 
previously. Qn account of the higher boiling 
point of the m-propyl nitrite and of some of the 
reaction products, it was necessary to heat the 
connecting tubing to a higher temperature than 
before. Stopcocks leading to the reaction vesse! 
were therefore replaced by mercury seals, and 
the whole of the connecting tubing was heated 
during experiments to 105°C. Temperature 
measurement and control was as before. It was 
possible to keep the temperature constant to 
within 0.25°C. 

Propyl nitrite was prepared in an analogous 
way to the previous nitrites, by the addition of 
a solution of nitrous acid to n-propyl alcohol. 
The liquid layer of nitrite was separated and 
purified by treatment with sodium carbonate 
solution, followed by several washings with 
ferrous sulphate solution to remove nitrogen 
oxides. It was then dried by permitting it to 
stand in contact with anhydrous sodium carbon- 
ate. The nitrite was then fractionally distilled 
and condensed twice. During the course of the 
investigation the n-propyl nitrite was stored in a 
cooled bulb. 
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DECOMPOSITION OF GASEOUS ALKYL NITRITES 


THE PRODUCTS OF THE REACTION 


Since the form of the reaction curves and total 
pressure increase at completion did not vary 
with temperature or pressure, it follows that 
there cannot be much variation in the products 
of reaction with varying conditions. 

Gas analysis was carried out on the products 
for nitric oxide as this was shown to be the only 
gaseous product when the reaction bulb was 
cooled down to —80°C and the gas sample 
pumped off for analysis. Duplicate analysis gave 
a value of 97 percent average for NO gas present. 
This result is in agreement with the overall 
reaction given below that is advanced for the 
decomposition : 


This equation requires that one molecule of 
NO be formed from one molecule of propyl 
nitrite. To confirm this, products were cooled to 
—80°C at which temperature the noncondensed 
NO gas registered the only pressure while the 
other products remained condensed. Instead of 
finding the NO pressure to be equivalent to that 
of the initial propyl nitrite, it was found to be 
invariably two-thirds. This low value is probably 
due to the formation of a complex between part 
of the NO and the condensed products at low 
temperatures. It may be stated here that a 
similar behavior was encountered in making runs 
where it was observed, with increasing initial 
pressure of propyl nitrite, there was a relative 
decrease in reaction rate. This difficulty was 
immediately overcome by raising the tempera- 
ture of the cooler connecting tubing to 105°C 
when all runs at one temperature gave constant 
rates independent of the initial pressure, as will 
be shown later. 

Thus it may be concluded from these facts 
and with a doubling of the pressure on comple- 
tion of reaction that the decomposition of 
n-propyl nitrite is analogous to the previous 
nitrites investigated where the initial break is of 


the O— NO bond. 
THE PrEssURE CHANGE ACCOMPANYING THE 
REACTION 


The mechanism given above would lead to a 
pressure increase of 100 percent at completion. 


TABLE I. Maximum pressure increases. 








MAXIMUM PRESSURE 
INCREASE (%) 


TEMPERATURE INITIAL PRESSURE 
(°C) (cm) 


189.9 
189.9 
189.9 
199.9 





26.68 
18.34 
14.68 
26.52 


98.4 
100.0 
99.2 
99.2 


Mean 99.2 








The pressure increases obtained were almost 
exactly 100 percent. Some typical values are 
given in Table I. 

We have therefore assumed that the pressure 
increase is a direct measure of the reaction, and 
that Ty2.5, T2; and Tso (the times for 12.5, 25, 
and 50 percent pressure increase) correspond to 
12.5, 25 and 50 percent reaction, respectively. 
On this basis the ratio 759/72; should be equal 
to 2.44 if the reaction is first order. Actually 
the mean value of the ratio from all runs was 
2.53. 


THE RATE OF REACTION 


Complete data for two typical runs are given 
in Table II. The first-order constants obtained 


TABLE II. Data for typical runs. 


Temperature = 190°C 


Temperature = 210°C 
Initial pressure = 14.68 cm 


Initial pressure = 18.38 cm 








TIME AP 7, 


TIME AP X 5 
(min.) (cm) REACTION K 


(min.) (cm) 


ar 
4/0 
REACTION K 





0.0257 
0.0272 
0.0273 
0.0257 
0.0237 
0.0226 
0.0227 
0.0227 
0.0223 
0.0227 
0.0228 
0.0227 
0.0223 
0.0217 
0.0216 
0.0214 
0.0216 
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are slightly high at the start, and then become 
very satisfactorily constant. 
The summarized data for all runs are given in 


Table III. 


THE EFFECT OF PRESSURE ON THE RATE 


The constancy of the values of the fractional 
times with varying pressure shows that the 
reaction is of the first order. There is no indica- 
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tion of any falling off in rate at the lower pres- 

sures used. The investigation of the nitrite 

decompositions at low pressures is in progress 

and the results will be reported later. 


THE EFFECT OF SURFACE 


Experiments were made in an empty reaction 
vessel, and also in a vessel packed with short 
lengths of Pyrex tubing, so as to give an increase 
of about 9-fold in the surface/volume ratio. It 
will be seen from Table III that there was very 


TABLE III. The rate of reaction. 








(em) Ti2.5 (see.) T2gec.) Tao(see.)| (°C) (em) Ti2.5(see.) T25(see.) T0(sec.) 





Empty bulb 
132.0 306.0 
148.8 
144.0 
150.0 
131.1 
51.0 
58.8 
54.0 
69.0 
67.8 
ea 60.0 
Packed bul 
9.32 265.2 
11.26 265.8 
18.75 


319.8 
Mean 283.6 
13.05 


122.4 
14.92 


121.2 
Mean 


121.8 
16.70 57.6 


Empty bulb 
7 1998 4740 
864 
936 


19.24 
26.52 
27.65 
31.24 
Mean 


199.9 


b 

666 
660 
774 
700 
297.6 
301.2 
299.4 
136.8 











little difference in rate in the two vessels. There 
is therefore no doubt that in the empty vessel 
there are no complications due to heterogeneity. 


THE TEMPERATURE COEFFICIENT 


In Fig. 1 the logarithms of the fractional times 
are plotted against 1/7. The data used are the 
mean values given in Table III. From the slopes 
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of the lines we obtain for the heat activation: 


37,550 calories per gram molecule 
37,550 
37,850 


from T12.5 
T 25 
Tso 





Mean 37,650 
Hence the rate of reaction may be given by 
T25 = 1.05 X 10-e3765/ FT sec, 
or, since KT; =0.288, 
K =2.75 KX 10'4e—37650/2T sec —1, 


DISCUSSION 


In the past attempts have been made to 
compare the unimolecular decompositions of 
series of related compounds. In most cases, 
however, such attempts have not been very 
successful on account of complications, such as 
side reactions, which make the decompositions 
of the different members of a series not strictly 
comparable. In addition, in most cases previously 
investigated two or more chemical bonds are 
being broken in the overall reaction process, and 
this makes the actual elementary reaction very 
uncertain. 

In the present case, there seems to be no 
doubt that all the nitrite decompositions so far 
investigated are directly comparable, and that 
in each case the primary reaction involves the 
rupture of only one chemical bond. The nitrites 
therefore furnish an ideal series for purposes of 
comparison. Obviously such a comparison can 
be more accurately made when measurement at 
low pressures have been completed, and when 
data are available for higher members of the 
series. It seems, however, to be worth while to 
make a qualitative comparison of the results so 
far obtained. 

The rates of the three decompositions are 
given by an equation of the form 


K=Ae-#!®" viz: 

K =1.84 X 10'3¢—36400/2T soc—!, 
K =1.39 X 10'4¢@-37700/RT csoc-!, 
K =2.75 XK 10'4e—37650/2T sec. 


Methyl Nitrite 
Ethyl Nitrite 
n-Propyl Nitrite 


The chief difficulty in making comparisons is 
that the exponential and nonexponential factors 
in the rate expression are interdependent. As a 
result experimental errors cause fluctuations in 
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the two which are much greater than the error 
in determining the absolute rate. It is therefore 
advisable perhaps, first of all to compare the 
velocity constants at some definite temperature. 
At 189.9°C we have 


Methyl 
Ethyl 
Propyl 


K=0.97 X10~ sec.“!, 
K=1.89X10~ sec.—!, 
K=3.95 X10 sec.-!. 


Whence we have 


Kure >:Kr :Kp,=1 Ps 1.95 : 4.07. 
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In other words the rate is almost exactly doubled 
on going from one member to the next. This is 
also equivalent to stating that, if we assume that 
the variations in E are due to experimental error, 
then A is approximately doubled as we go from 
one member of the series to the next. It seems, 
however, very unlikely that E should remain 
constant while A varies in this way, since the 
variation in the collision number will be com- 


' paratively small. In any case, however, there is 


no doubt that E changes only very slightly as 
we go up the series. 
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The Constancy of Two Dipole Moments in the Vapor State and Their Apparent 
Variation in Solution 


C. P. SMytu AND K. B. MCALPINE, Frick Chemical Laboratory, Princeton University 
(Received March 23, 1935) 


The validity of the solution method for the determination 
of dipole moment is considered with particular reference 
to the possible variation of moment with temperature. 
The moments of heptyl bromide and butyl. chloride in the 
vapor state are found to be 2.07+0.05 and 2.04+0.01, 
respectively, which show the absence of appreciable in- 
crease in moment in the normal alkyl chlorides and 


HESE measurements in the vapor state 

were made for the sake of careful com- 
parison with previous measurements in solution 
upon two substances which apparently showed 
abnormally and improbably high values of the 
atomic polarization. The comparison takes on 
added interest now because of recent work on 
the effect of solvent upon the apparent value of a 
dipole moment determined in solution. 

Because of the rise of the apparent value of 
the dipole moment determined in solution with 
rising temperature, Jenkins' has stated that 
measurements in solution cannot be used for 
studying the variation of the real moment with 
temperature. Obviously, a small variation of the 
apparent moment in solution with rising temper- 
ature is not to be taken as evidence of variation 
in the real moment of the molecule, nor has it 
been generally so regarded. Many apparent 


'H. O. Jenkins, Trans. Faraday Soc. 30, 739 (1934). 


bromides with increase in the length of the carbon chain 
beyond ethyl. Consideration of these values together with 
data for the alkyl halides in solution leads to the following 
estimated moment values for the vapor state: n-propy] 
bromide and n-butyl bromide, 2.04+0.03; n-propyl iodide 
and n-butyl iodide, 1.93+0.04. 


small rises in moment from 25° to 50° in solution? 
have not been treated as significant. An ab- 
normally large rise, however, in the case of 
diethyl succinate led to further investigation® 
which established the reality of the effect along 
with the rise of moment with temperature in 
ethylene chloride and chlorobromide. At the 
same time, an apparent smaller rise for ethyl 
acetate was disregarded. 

Obviously, as the real variation of moment 
with temperature depends upon the potential 
energy between the movable parts containing 
the dipoles, the presence of a solvent around the 
molecule may be expected to affect the effective 
value of the moment as well as its variation. 
Thus, the moments found for ethylene chloride, 


2C. P. Smyth and W. S. Walls, J. Am. Chem. Soc. 53, 
527, 2115 (1931); ibid. 54, 1854, 3230 (1932). 

8C. P. Smyth, R. W. Dornte and E. B. Wilson, Jr., 
J. Am. Chem. Soc. 53, 4242 (1931). 
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chlorobromide and bromide*:‘* in heptane solu- 
tion are somewhat higher than those found for 
the vapors of these substances at the same 
temperatures® and show a 50 percent greater 
increase with temperature. The moments in 
benzene solution are much higher and almost 
independent of temperature, presumably, be- 
cause of some specific effect of the solvent.® 
Ethylene chloride in solution in carbon disulfide, 


carbon tetrachloride and chloroform has moment ° 


values very slightly lower than in heptane® and 
with temperature coefficients like that of the 
vapor. Because of the tendency toward variation 
with temperature of apparent moment values in 
solution, small differences in these temperature 
coefficients cannot be regarded as significant. 
Since no relation has been found between ‘the 
apparent moment of ethylene chloride in solu- 
tion®: 7 and the dielectric constant of the solvent, 
the moment being higher in solution than in the 
vapor state instead of lower as usually found, it 
is impossible to apply to these molecules with 
movable dipoles any of the empirical relations 
which have been developed between apparent 
moment and dielectric constant of the solvent. 
In the case of such molecules, it is obvious that 
the variation of polarization with temperature 
cannot be used to calculate the atomic polariza- 
tion. Such values calculated by Smyth and 
Kamerling* for four polymethylene bromides 
containing two movable dipoles were enclosed 
in parentheses and taken as evidence of variation 
of the moment with temperature which made the 
calculated atomic polarization impossibly high. 
These values should not, therefore, be cited as 
evidence of the incorrectness of the temperature- 
solution method of obtaining moments.! 


PURIFICATION OF MATERIALS 


n-Heptyl bromide 


Material from the Eastman Kodak Company 
was purified as in earlier work;® mp?® 1.45058. 





4C, P. Smyth and S. E. Kamerling, J. Am. Chem. Soc. 
53, 2988 (1931). 
: aaa T. Zahn, Phys. Rev. 38, 521 (1931); ibid. 40, 291 
(1 R 

6 A, E. Stearn and C. P. Smyth, J. Am. Chem. Soc. 56, 
1667 (1934). 

7H. Miiller, Physik. Zeits. 34, 689 (1933). 

8 C. P. Smyth and H. E. Rogers, J. Am. Chem. Soc. 52, 
2227 (1930). 
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Two samples were used in the dielectric constant 
measurements, identical results being obtained. 


n-Butyl chloride 


Material from the Eastman Kodak Company 
was refluxed gently with concentrated sulfuric 
acid, washed several times with water, dried 
over two portions of calcium chloride and 
distilled ; b.p. 78.4-78.6° ; np”® 1.40230. Repeated 
treatment did not change the refractive index. 


EXPERIMENTAL RESULTS 


The dielectric constants « were determined 
with the apparatus previously described,® hepty] 
bromide being measured in a copper cell and 
butyl chloride in a gold-palladium cell. The 
polarization P was calculated by means of the 
equation P=[(e—1)/(e+2)]V, in which V is 
the molar volume. For heptyl bromide, the 
polarization at each temperature was taken as 
the average of from three to eight measurements 
at pressures between 35 and 162 mm, the error 
caused by deviation from the ideal gas laws at 
these pressures being small. For butyl chloride, 
five to eight determinations at each temperature 
were made over a range of pressure extending to 
higher values and extrapolated to zero pressure 
to obtain a value Py which should be free from 
error due to gas law deviations. For hepty! 
bromide the temperature range of the measure- 
ments was so limited by the tendency toward 
thermal decomposition at higher temperatures 
that the dipole moment yw was calculated at 
each temperature by means of the equation 
w=0.01273[(P—MRp)T |}, in which MRp is the 
molar refraction for the D sodium line. The 
more accurate equation P=a+b/T was applied 
to the butyl chloride results and the moment was 
calculated as »=0.01273X10-%5!, while the 
atomic polarization was obtained as P,4=a—Pz, 
the difference between a, the total induced 
polarization, and Pz, the molar refraction extra- 
polated to infinite wavelength. 


DISCUSSION OF RESULTS 


The two values previously found for the 
moment of heptyl bromide, 1.83 in heptane 


°K. B. McAlpine and C. P. Smyth, J. Am. Chem. Soc. 
55, 453 (1933). 










n= — = — > rr 


—_—[—S >|, —- —§ Pp -&- — TH —-—- —-— -— —- |. FF KRY 


—_—_ re fF -—> =F . 


1wo we 


— 





DIPOLE MOMENTS 


solution® and 1.85 in benzene solution,'® were in 
excellent agreement, but the value obtained in 
heptane solution was calculated from the varia- 
tion of the polarization with temperature. This 
method gave the high value 15.5 for the atomic 
polarization, values of this magnitude having 
since been shown to be incorrect," and a corre- 
spondingly low value for the moment. When 
the moment is calculated from the value obtained 
by Smyth and Rogers for the polarization at 
each temperature by means of the equation 
u=0.01273[ (P.,— MRp)T }}, the values rise from 
1.97 at 223°K to 2.02 at 303°K. One may apply 
to these data the empirical equation of Miiller,’ 
Pyuges= Peoi/{1 —0.075(e— 1)? ], where Pu 
=P.,—MRp and « is the dielectric constant of 
the solvent. The moment values obtained from 
these values of Py calculated for the gaseous 
condition increase from 2.05 at 223°K to 2.09 
at 303°K. If 3.5 is added to the value 42.3 used 
for MRp in order to account for the atomic 
polarization, the moment calculated for the gas 
from the liquid measurements increases from 
2.02 to 2.04 over this range of temperature. 
The use of this value Pg +P,4=45.8 to calculate 
the moment from the gas measurements gives a 
mean value 2.09. Since this is higher than the 
mean value 2.03 calculated in the same way from 
the empirically corrected solution values and 
since such deviations from the gas laws as may 
have occurred would tend to make the polariza- 
tion and, hence, the moment values high, it 
seems reasonable to assign a value 2.07+0.05 to 
the moment of heptyl bromide. The difference 
between this value and the accurate value 2.02 
found for ethyl bromide vapor” is equal to the 
probable error in the determination. 

The accurate values for butyl chloride in 
Table I show, as in the case of heptyl bromide, 
that the atomic polarization value 13.6 obtained 
from the variation of polarization in heptane 
solution with temperature® was high and the 
moment value 1.88 correspondingly low. Calcu- 
lation of the moment from these solution polari- 
zations by subtraction of the molar refraction 
gives values rising from 1.98 at 183°K to 2.07 


'0 J. Errera and M. L. Sherrill, J. Am. Chem. Soc. 52, 
1993 (1930). 

'C. P. Smyth, J. Chem. Phys. 1, 247 (1933). 

"C. P. Smyth and K. B. McAlpine, J. Chem. Phys. 2, 
499 (1934), 


IN VAPORS AND SOLUTIONS 


TABLE I. Polarizations and dipole moments. 








n-Butyl chloride 


n-Heptyl bromide 
i, = 4 T, °K Po 


ux 10'8 





372.8 
373.5 
377.0 
388.5 
393.0 
395.5 
414.8 
415.8 
423.3 
433.8 


2.13 
2.17 
2.13 
2.14 


314.7 
337.4 
379.1 
405.5 
nae 428.2 88.15 
2.17 479.6 82.0; 
2.16 a=28.5, P4=3.6, »=2.04 X 10718 
2.15 
, 2.16 

108. 2.15 
Mean value =2.15 


109.9 
104.5 
96.4 
91.7 
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at 343°K. These values are close to those found 
by Parts," 1.97 at 283°K in benzene solution and 
the approximate value 2.00 calculated from a 
variation of polarization with temperature over 
a range of only 40°. Calculation of the polariza- 
tions in heptane to the vapor state by means of 
the Miiller equation leads to moment values 
rising from 2.07 at 183°K to 2.13 at 343°K. If, 
however, the polarization values thus calculated 
to the vapor state are diminished by 3.0, the 
difference between the value of a in Table I and 
MRp, the moment variation is reduced to a rise 
from 2.05 to 2.09 over this temperature range. 
These values are close to the accurate value 
2.04 in Table I. 

The identity of the moment of ”-butyl chloride 
with the accurate value 2.04 found by Sanger" 
for n-propyl chloride and the closeness of both 
values to the most accurate found for ethyl 
chloride, 2.05 by Fuchs'® and 2.02 by Sanger 
shows the absence of detectable increase in 
moment with increase in the length of the chain 
of the alkyl chloride beyond two carbons. The 
variation with temperature of the polarization of 
butyl bromide in solution in heptane® gave a 
high atomic polarization 19.3 and a low moment 
1.81. Calculation of the moment by subtracting 
the molecular refraction MRp gives 1.97 at 
183°K, identical with the value found by Parts 
at 283°K in benzene solution, and 2.08 at 343°K. 
Subtraction of an additional 3.5 to take care 
approximately of the atomic polarization and 
application of the Miiller equation to these data 


13 A. Parts, Zeits. f. physik. Chemie B7, 327 (1930). 

4 R, Sanger, Helv. Phys. Acta 3, 161 (1930); Physik. 
Zeits. 32, 20 (1931). 

15Q,. Fuchs, Zeits. f. Physik 63, 824 (1930). 
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give 2.02 at 183°K and 2.09 at 343°K. The 
mean of these two values 2.06 is indistinguishable 
from the mean 2.07 found in similar fashion for 
butyl chloride. With the moment of ethyl 
bromide vapor 2.02," it is probable that the 
moments of propyl and butyl bromide vapors 
differ by no more than 0.03 from 2.04, which 
differs by no more than the probable error from 
the value 2.07 here determined for heptyl 
bromide. The variation with temperature of the 
polarization of butyl iodide in solution® gave a 
very high atomic polarization 21.9 and a corre- 
spondingly low moment value 1.59. Calculation 
of the moment from these data by use of the 
molecular refraction gives values rising from 1.80 
at 193°K to 1.87 at 293°K. In view of the 
acturate value 1.90 for the moment of ethyl 
iodide vapor,” it is probable that the moments 
of n-propyl and of n-butyl iodide vapor differ 
from 1.93 by no more than 0.04. 

It is evident from these calculations of ap- 
parent moment values at different temperatures 
in solution as well as from a tendency toward 
increase between values at 25° and 50° previously 
published? that, in many cases, the apparent 
value of the moment in solution increases with 
rising temperature. The change of polarization 
with temperature thus involved affects the slope 
of the P—1/T curve so as to give too large a 
value of a and a correspondingly small value of b. 
Although the method of plotting polarizations 
in solution against 1/7 has given satisfactory 
values for the moments in solution of chloro- 
benzene, chloroform! and f¢ert.-butyl chloride,!” 
it is evidently wholly unsatisfactory for the 
n-alkyl halides and cannot, therefore, be gener- 
ally employed with safety. 

In view of the dependence of the apparent 
polarization and moment of the solute upon the 
dielectric constant of the solvent found by 


16C, P. Smyth and S. O. Morgan, J. Am. Chem. Soc. 
50, 1547 (1928). 

17C, P. Smyth and R. W. Dornte, J. Am. Chem. Soc. 
53, 545 (1931). 
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Miiller,’: '® it follows that the apparent value 
of a moment in solution should rise as the 
dielectric constant of the solvent decreases with 
rising temperature. Jenkins! has been able to 
account for increase in the apparent moment of 
nitrobenzene in solution in dekalin with rising 
temperature entirely by the effect of the decrease 
in the dielectric constant of the dekalin. How- 
ever, the results of the calculations given in this 
paper show that in the case of heptyl bromide 
and butyl chloride, the effect of change in the 
dielectric constant of the solvent as calculated 
directly by Miiller’s equation accounts for only 
about one-third of the apparent variation. When 
allowance is made also for atomic polarization, 
the apparent variation is reduced to about 
one-half. The qualitative explanation previously 
given!! for the incorrectness of the abnormally 
high atomic polarization values involves much 
the same fundamental idea as do these empirical 
equations for the effect of solvent. It may be 
said that the apparent moment found for a 
molecule in solution is influenced by the internal 
field acting on the molecule. Although the effect 
of this internal field, which decreases with rising 
temperature, is greatly reduced by extrapolation 
to infinite dilution and, sometimes, wholly elimi- 
nated by the further use of an empirical equation 
involving the dielectric constant of the solvent, 
a small error may yet remain in the moment 
values obtained from solution. As a matter of 
fact, this was shown by Miiller’s discovery that 
there were other substances in addition to those 
containing movable dipoles to which his equation 
did not apply. The failure to eliminate entirely 
this effect of the internal field may evidently 
have only a small influence upon the moment 
value determined at one temperature, but may 
seriously impair the accuracy of values deter- 
mined from the variation of polarization with 
temperature. 


18 See also, H. O. Jenkins, Nature 133, 106 (1934) and 
S. Sugden, ibid. 133, 415 (1934). 
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No satisfactory crystal structure determination has 
been made of any of the forms of phosphorus. Five samples 
of black phosphorus gave identical powder patterns 
different from those reported by Linck and Jung, whose 
rhombohedral structure gives improbable interatomic 
distances and coordination. A sample which had been 
prepared at room temperature and 35,000 atmospheres 
gave the diffuse rings of an ‘‘amorphous”’ x-ray pattern. 
Another sample, prepared at 300°C and 8000 atmospheres, 
is a new form of phosphorus, having the same density as 
black. Atomic distribution curves of crystalline and 
“amorphous” black and red phosphorus were obtained by 
the method of Fourier analysis. All four gave practically 


the same curve, showing three neighbors at 2.28A and 
about twelve at 3.6A. With the aid of the atomic dis- 
tribution curves and the fact that certain lines showed a 
preferred orientation, it was possible to determine the 
structure of black phosphorus. It consists of double layers; 
the cell is side-centered orthorhombic with a=3.31A, 
b=4.38A and c=10.50A. There are eight atoms in the cell, 
giving a calculated density of 2.69 compared with 2.70 
observed. The space group is V,'!*—Bmab. The black 
phosphorus structure allows a closer packing than the 
arsenic structure, where the atomic environment is similar, 
and so is favored by the high pressure. 





INTRODUCTION 


HOSPHORUS crystallizes in a number of 

allotropic forms, three of which are easily 
distinguished because of their widely differing 
physical and chemical properties. These three 
forms are the red, white and black varieties. 
Their peculiar conditions of formation and trans- 
formation make their study of interest in the 
theory of the dynamics of crystalline changes. 
The data on phosphorus are confusing and in- 
complete.! None of the crystal structures have 
been satisfactorily worked out by x-ray methods; 
the various investigators do not even agree on 
the positions of the lines in the powder patterns. 


PROPERTIES OF THE ALLOTROPIC FORMS OF 
PHOSPHORUS 


When heated, all three forms sublime to a 


vapor composed of P, molecules, which yields. 


white phosphorus when condensed. However, at 
lower temperatures (100°C in a vacuum), the 
vapor of red phosphorus condenses to give red 
phosphorus; but the vapor density is so low that 
no molecular weight determination has been 
made. White phosphorus is easily changed to 


* National Research Fellow. 
‘For a discussion of physical and chemical properties of 
phosphorus see Ephraim. Anorganische Chemie, Steinkopft, 
Dresden and Leipzig, 1934, pp. 87-92. 


red by the action of heat, light or x-rays. Under 
a pressure of 12,000 atmospheres and at a 
temperature of 200°C, it changes to black 
phosphorus.? Red phosphorus does not change 
under these conditions.t Black phosphorus 


changes to red when heated to 550°C. All of 


these changes are irreversible. 
The physical differences between these forms 
are striking as shown in Table I. 


TABLE I. Physical properties of the forms of phosphorus. 








M.P. B.P. INFLAMMABILITY SOLUBILITY 


44.1 280 


Cotor DENSITY 


White 1.83 





Very high Soluble in number 
of organic sol- 
vents with mol. 
wt. Py 
Insoluble 
Insoluble 


High 
Low 








These facts seem to indicate that white 
phosphorus is composed of more or less loosely 
bound P, molecules which readily break up to 
form the more stable structures of red and black 
phosphorus, which consist of larger molecules or 
else are continuous, nonmolecular structures. 

White phosphorus can be obtained in beautiful, . 
transparent crystals of the cubic system by 
crystallization from solution or by slow sublima- 


2 P. W. Bridgman, J. Am. Chem. Soc. 36, 1344 (1914). 

+ Professor Bridgman has recently succeeded in changing 
red phosphorus to black by the application of a shearing 
stress at very high pressures. 
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tion. Jung’ reports that x-ray analysis is im- 
possible since the x-rays transform the white to 
the red form. However, Natta and Passerini‘ 
announced that by cooling the sample to —35°C 
and by using iron radiation, they were able to 
obtain 22 lines of a cubic pattern with a2=7.17A. 
But they failed to publish their data or the 
crystal structure. 

Red phosphorus in its commercial form gives 
diffuse rings in its x-ray pattern and is supposed 
to be amorphous. After suitable heat treatment 
it gives a normal line pattern. Red phosphorus is 
usually considered to be either a finely divided 
form of violet phosphorus or a solution of white 
in violet because of similarity in physical 
properties. Frost,5 however, has shown by x-ray 
evidence that they have different crystalline 
forms. He concluded from his powder patterns 
that both structures had a low symmetry. 
Violet phosphorus can be obtained in small, 
transparent, monoclinic crystals by crystalliza- 
tion from molten lead. Olshausen® and Jung* 
report patterns which do not agree with either 
of the patterns of Frost, nor do they agree with 
one another. Olshausen used red phosphorus and 
Jung, violet. 

Black phosphorus resembles graphite when 
pulverized. Linck and Jung’ assign it a rhombo- 
hedral cell on the basis of their x-ray pictures. 
However, the coordination of atoms in their cell 
(one nearest neighbor at 2.46A and six at 3.01A) 
does not agree with the accepted atomic radius 
of 1.10A, nor does it give the three or five nearest 
neighbors to be expected from a _ covalent 
structure. 


EXPERIMENTAL WORK 


We were fortunate in securing seven samples ot 
black phosphorus from Professor Bridgman at 
Harvard. Five of these were prepared at the 
usual temperature of 200°C and pressures of 
12,000 atmospheres or more. These samples 
_varied in age from a few months to twenty 
years. All gave identical powder patterns which 


3 Jung, Centralblatt f. Min. u. Geol. 1926, 107-114. 

4 Natta and Passerini, Nature 125, 707 (1930). 

5 Frost, J. Russ. Phys. Chem. Soc. 62, 2235 (1930). 

6 Olshausen, Zeits. f. Krist. 61, 463 (1925). 

7 Linck and Jung, Zeits. f. anorg. Chemie 147, 288 (1925). 
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were not the same as those reported by Linck 
and Jung.’ Heat treatments of 100 hours at 
450°C and two hours at 550°C failed to change 
the patterns in any way except to make the lines 
somewhat sharper. Linck and Jung had used 
samples from the same source and the physical 
properties they describe are the same as those 
we observed. We are at a loss to account for the 
discrepancy in the patterns. 

The other two samples had been made under 
different conditions. One was made at room 
temperature and 35,000 atmospheres.* Its phys- 
ical properties were similar to those of the other 
samples except that it had a density of 2.6 
instead of 2.7. Its powder pattern showed diffuse 
rings instead of sharp lines. The other had been 
made by Dr. F. Birch at 300°C and 8000 
atmospheres and gave a powder entirely different 
in appearance from the other samples. It had a 
deep red color and no flakiness. Its x-ray pattern 
was different from any of the other forms of 
phosphorus. It is a new form of phosphorus 
which we have not found time to investigate 
further. 

Samples of commercial, “amorphous” red 
phosphorus were obtained from the stock room 
and showed diffuse rings on the powder pattern. 
We also obtained a sample of red phosphorus 
from Professor Bridgman which had withstood a 
pressure of 12,000 atmospheres without change, 
and prepared other samples by heating black 
phosphorus to 550°C in a sealed quartz tube. 
All of the samples except, of course, commercial 
red phosphorus, gave the same set of x-ray lines 
which checked that of Frost® for red phosphorus. 

We have applied the method of Fourier 
analysis to the x-ray powder patterns of black 
and red phosphorus. This method, which was 
developed by Zernike and Prins® and by Debye 
and Menke,!® yields information about the 
distances between atoms, but nothing about 
their arrangement in the crystalline lattice. It 
has been used as an aid in crystal structure de- 
terminations by Warren and Gingrich" and by 
Warren.” 


8 P. W. Bridgman, Phys. Rev. 45, 844 (1934). 

® Zernike and Prins, Zeits. f. Physik 41, 184 (1927). 

10 Debye and Menke, Erg. d. Techn. Réntgenkunde II. 
1 Warren and Gingrich, Phys. Rev. 46, 368 (1934). 

12 Warren, J. Chem. Phys. 2, 551 (1934). 





CRYSTAL STRUCTURE OF 


ATOMIC DISTRIBUTION CURVES 


Assuming that each phosphorus atom has the 
same environment, the density of atoms at a 
distance r from a given atom can be obtained 
from the following formula :" 


Arr p(r) =4 rr? po+ (27/7) fo°si(s) sin (rs)ds, 


s=47n sin 0/X, 

po=average density of sample in atoms per unit 
volume, 

p(r) =average density of atoms at a distance 7 from 
a given atom, 
i(s)=(I/N—f?)/f*, 

J =experimental intensity corrected for absorp- 
tion, incoherent radiation and polarization, 
including polarization due to the mono- 
chromatizing crystal, 

f?=intensity of independent scattering per atom, 

N=number of scattering atoms. 


where 


Consequently 4zr’p(r)dr represents the number 
of atoms between r and r+dr. 

The analysis was carried out in the usual 
way! except for the determination of the 
constant N, which was evaluated from the 
following formula :* 


N=f0'Isds/fo'f?sds. 


Molybdenum Kea radiation was used, mono- 
chromatized by a sodium chloride crystal. The 
camera had a radius of 47 mm and the exposures 


* At large angles the experimental error in the determi- 
nation of intensity is high, partly due to the small in- 
tensities involved and partly due to the large uncertain 
correction for incoherent intensity. Hence the assumption 
that J= Nf? at some arbitrary large angle may give a 
result which is far from correct. A better way is to choose 
N so that the deviations of the experimental curve from 
independent scattering more or less balance one another 
at large angles. A method which depends only a little on 
the inaccurate measurements at large angles is the fol- 
lowing. From the formula :" 


T=NP {tS 42rLo(r) —pol(sin (sr)/sr)dr} 


we easily obtain 


Jf tsds -N f Spsds 
=42N{" {*p sin (sr)r[p(r) —po jdsdr = K. 


It can be shown that if S is chosen large enough, the value 
of K is small compared with that of the other two integrals, 
80 it is justifiable to assume K to be zero as a first approxi- 
mation. We found after doing this that its actual value 
was about eight percent of the other integrals, hence a 
second approximation was unnecessary. The error in the 
first approximation can be estimated in advance by inte- 
grating K over the region where p(r) is known to be zero. 
In the case of phosphorus this estimate was 11 percent. 
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Fic. 1. Atomic distribution in black phosphorus. The 
solid vertical lines indicate the distribution in the crystal 
structure found for black phosphorus. 
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Fic, 2. Atomic distribution in red phosphorus. 
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required several days. It was found that errors 
due to stray radiation or to incorrect determina- 
tion of N had a considerable effect on the values 
of density determined, but that the positions of 
the peaks remained remarkably constant. 

The results of the analysis are shown in Figs. 
1 and 2. The four samples analyzed; crystalline 
and ‘‘amorphous” black phosphorus and crystal- 
line and commercial red phosphorus gave nearly 
the same distribution curves. The supposedly 
amorphous samples each show four bands at the 
positions where the crystalline samples have 
their strongest lines. Probably they are simply 
extremely fine crystals. The area under the first 
peak is 3.3, indicating three neighbors at about 
2.28A. The second peak is at 3.6A and has an 
area indicating about twelve atoms at that 
distance. The red phosphorus distribution curves 
show peaks slightly farther out than those of 
black. 


THE CRYSTAL STRUCTURE OF BLACK 
PHOSPHORUS 


The interplanar spacings of black phosphorus 
were accurately determined from powder pat- 
terns made with filtered iron radiation from a 
Hadding tube. In one case iron and in another 
sodium chloride was mixed with the sample as a 
standardizing substance. The results usually 
varied from one another less than 0.005A. 
Supplementary photographs from a copper tube 
checked the intensities and approximate positions 
of the lines, and also yielded a few more faint 
lines. The camera used had a radius of 47 mm. 

The interplanar spacings found will not fit a 
cubic, tetragonal or hexagonal lattice. Ordinarily, 
it is hopeless to try to solve a structure of lower 
symmetry from powder photographs alone. 
Fortunately, the results of the Fourier analysis, 
together with the fact that one of the lines and 
its two higher orders showed preferred orienta- 
tion, made it possible to determine the structure. 
The physical properties of the powder and the 
preferred orientation of the lines indicated a 
layer structure with the layers 5.25A apart. 

The second order of the line showing preferred 
orientation is stronger than the first, while the 
third is weaker. A rough calculation shows that 
each layer must be double with a spacing of 
about 2.1A between halves. In satisfying this 
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condition and the distribution determined from 
the Fourier analysis, we are led to an ortho- 
rhombic cell of which we know the approximate 
dimensions. Trial shows that the c axis must be 
double the layer spacing and more accurate 
values for the a and b axes are determined from 
the first few lines. It soon becomes evident that 
the true lattice is side-centered orthorhombic 
with the following dimensions: 


a=3.31A; b=4.38A; c=10.50A. 


The unit cell contains eight atoms, which gives a 
calculated density of 2.69, compared with the 
observed value, 2.70. The space group is 


V;,:8— Bmab 


with the eight atoms in the following positions: 


ei il el 1 ° 1 1 . 
Ouv; 2, 2—U, U; 2, U, at; 0, 2—4U, a+; 
Oud; 2, 2+u, 0; 2, H, 2-0; 0, 3+u, 3-v. 


The best agreement of intensities is obtained for 
u=0.090 and v=0.098. 

Calculated and observed interplanar spacings 
and intensities of the first 42 possible lines of 
this space group are shown in Table II. The 


TABLE II. Calculated and observed interplanar spacings and 
intensities for black phosphorus. 








deale. Teale. Iobs.| hkl deale. dobs. 


Teale. Tobs. 





180 214 1.334 


1.33* 
131 1.325 — 
220 = 1.320 — 
008 = 1.313 
117 1.304 
222 : — 
034 ; 1.27* 
018 25 1.249 
133 248 — 
206 203 —-:1.200 
224 j 1.18* 
216 f 
127 “a 1.156 
135 f — 
028 ] — 
0361. — 
040 ~—-1.095 — 
232 
204 042 
125 J _ 119 
026 = 1. 1.364 j 3/311 


1.30* 
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032 
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1.065 











* Measured with less accuracy from copper radiation. 


agreement is excellent in nearly all cases. The few 
small deviations are due to preferred orientation, 
which cannot always be allowed for, tending to 
make planes near (00/) stronger and other planes 
weaker. Allowance for the temperature factor 
would make the calculated intensity for planes of 
smaller d, smaller. After the unit cell is once 
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discovered, the determination is rigorous. The 
atomic distribution is shown in Fig. 1 for the 
crystal as well as that determined by Fourier 
analysis. The agreement is very good, except 
that the Fourier analysis gives slightly too large 
distances. 


DISCUSSION OF RESULTS 


The unit cell of black phosphorus consists of 
two double layers, one of which is shown in 
Fig. 3. Each atom is bonded to three nearest 
neighbors at 2.18A, in agreement with the 
accepted atomic radius of 1.10A, and with the 
coordination to be expected from covalent bond- 
ing of phosphorus. Two of the bonds are in the 
plane of the layer, at 99° from one another; 
the third is between layer halves at 103° 30’ from 
both, making the average bond angle 102°. 
This agrees with the tendency of bond angles to 
be nearer tetrahedral for the lighter elements of 
the periodic group. Thus for the fifth period the 
trend is: P 102°, As 97°, Sb 96° and Bi 94°. 
The atomic environment is quite similar to that 
in the arsenic structure, which is assumed by the 
other members of the group, except nitrogen. 
The difference consists in the fact that in the 
latter structure all bonds are between halves 
of the double layer, while in the former only one- 
third of the bonds are. 

The question arises as to why phosphorus does 
not assume the arsenic structure. A good reason 
is that the black phosphorus structure is closer 
packed than the arsenic structure and so is 
favored by the high pressure under which black 
phosphorus is formed. If phosphorus assumed 
the arsenic structure, retaining the bond dis- 
tances and angles of black phosphorus and also 
the closest distance of approach between atoms 
in different layers, it would have a density of 
only 2.44 instead of 2.69. Arsenic at high 
pressures might be expected to take the black 
phosphorus structure since, by doing so, its 
density would be increased from 5.73 to 6.43. 


STRUCTURE OF 


BLACK PHOSPHORUS 
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F1G. 3. A double layer of the black phosphorus structure. 
‘““4+" and ‘‘—” indicate that the atom is, respectively, 
0.098 c above or 0.098 c below the plane of the paper. 


The fact that it has not been observed is not a 
conclusive argument against its possibility, since, 
as we saw, black phosphorus forms only under 
special conditions. 

In this connection it should be pointed out 
that arsenic is a good deal more metallic than 
black phosphorus, having a specific resistance of 
35 X 10-® ohm cm compared with 0.711 ohm cm 
for the latter.? The decisively covalent character 
of phosphorus is nicely shown by the fact that 
the distance between bonded atoms is only 
2.18A compared with 3.68A as the closest 
approach between atoms in different layers. In 
arsenic the bonds are less sharply differentiated ; 
three nearest neighbors are at 2.51A, with three 
more atoms at 3.15A. 

The Fourier analysis furnishes a valuable clue 
to the structure of red phosphorus. We know 
that each atom forms three covalent bonds in 
the normal manner. The fact that the second 
peak is at practically the same position as in 
black phosphorus suggests that the bond angles 
are nearly the same. 
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The Number of Canonical Structures of Each Degree of Excitation for an Unsaturated 
or Aromatic Hydrocarbon 


G. W. WHELAND, Gates Chemical Laboratory, California Institute of Technology 
(Received March 5, 1935) 


A simple procedure is outlined for obtaining the number of canonical structures of each 
degree of excitation for an arbitrary hydrocarbon molecule. The results are applied to the 
approximate calculation of the energies of a number of such systems. 





N the quantum-mechanical treatment of the 

unsaturated and aromatic hydrocarbons, and 
of the free radicals, it is frequently desirable to 
know the number of canonical structures of 
each degree of excitation.? These figures can be 
obtained, of course, by writing down all ‘the 
structures, and by then determining the number 
of each type by actual count. This becomes 
quite impracticable, however, with complicated 
molecules,*? and some simpler procedure is neces- 
sary. 


DESCRIPTION OF THE METHOD 


The orbitals are first arranged formally in a 
circle in the usual manner.‘ There is some 
arbitrariness in this step, since in general the 
final result will depend upon the way in which 


1L. Pauling and G. W. Wheland, J. Chem. Phys. 1, 
362 (1933). L. Pauling and J. Sherman, ibid. 1, 679 (1933). 
G. W. Wheland, ibid. 2, 474 (1934); 3, 230 (1935). 

2 The significance of the terms used here is as follows 
(cf. reference 1): In the method which is used by Pauling 
and his co-workers, and to which the conclusions reached 
in this paper are applicable, only the interactions between 
the [p ], orbitals are considered. These latter are arranged 
formally in a circle or polygon, and all structures are 
drawn in which each bond lies wholly within the figure, 
and in which no two bonds intersect. The structures thus 
obtained form the ‘‘canonical” set, and are the only ones 
which need be considered (cf. reference 4). The various 
bonds are of two types: ‘‘effective’’ bonds between orbitals 
that are actually adjacent in the molecule, and ‘‘ineffective”’ 
bonds between orbitals that are not adjacent. To the 
degree of approximation of the present treatment, the 
latter do not contribute to the stability of the molecule 
and are introduced solely in consequence of the mathe- 
matical formalism of the method. A structure in which 
all bonds are effective is called ‘‘unexcited’”’; a structure 
in which one bond is ineffective is called ‘‘first excited,” 
and so on. 

3 This is seen from the fact that the number of canonical 
structures for a system of 2N orbitals, each occupied by a 
single electron, is given by the expression (2N)!/N!(N+1)!, 
which increases very rapidly with N. Thus, with 2N=14, 
there are already 429 structures. 

4G. Rumer, Nachr. d. Ges. d. Wiss. zu Géttingen, 
M. P. Klasse, p. 337, 1932. L. Pauling, J. Chem. Phys. 1, 
280 (1933). 


the initial ordering is carried out. Thus, if the six 
orbitals in benzene are arranged in the sequence 
1, 2, 3, 4, 5, 6, in which they come in the molecule 
itself, there are 2 unexcited and 3 first excited 
structures. On the other hand, if they are 
arranged in the sequence 1, 3, 5, 2, 4, 6, there 
are 0 unexcited, 1 first excited, 3 second excited, 
and 1 third excited structures. This causes no 
great difficulty, however, since there is usually 
one single order which is more logical than any 
of the others (for benzene, for example, the first 
of the two given above), and, in any case, the 
problem is uniquely determined when once the 
order has been established. 

The general procedure will be to associate with 
each molecule, or grouping of atoms, a poly- 
nomial of the form 


Rot+kiz+ho2?+hy2° +--+ +hkn2", 


where k; is the number of structures of the ith 
degree of excitation, and 2 is simply a parameter. 
Obviously, finding the polynomial which corre- 
sponds to a given molecule is completely equiva- 
lent to solving the problem as originally formu- 
lated, since the coefficients k; are themselves the 
desired quantities. 


I. THE MOLECULE CONSISTS OF A STRAIGHT 
CHAIN 


The orbitals are arranged in the circle in the 
same order in which they come in the chain. 
The polynomial Cy, associated with the chain of 
2N atoms, is given by the recursion formula 


N 
Cy =Cy-1+22,C;-1Cw-,;, 


j=2 


which is derived in the following manner: 


356 
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We first number the orbitals in order, from 1 
at one end of the chain to 2N at the other end. 
Now we consider the canonical structures which 
contain a bond between orbital 1 and some other 
definite orbital, say 27. This bond divides the 
molecule into three parts: (1) the orbitals 1 and 
2j involved in the bond; (2) the orbitals 2, 3, 

- 27—1 lying on one side of the bond; and 
(3) the orbitals 27+1, 27+2, --- 2N lying on the 
other side of the bond. (In the special cases in 
which j=1, or N, the second, or the third, part 
will be missing. We then simply treat the missing 
part formally as a chain of 0 orbitals, and 
associate with it the polynomial C)=1.) If these 
three parts are considered separately, they can 
be represented by the polynomials J;, Jz and J3, 
respectively. Then J;=1 if 7=1 (i.e., if the bond 
between orbitals 1 and 2) is effective) and J,;=2z 
otherwise (i.e., if the bond is_ ineffective). 
Similarly Je=Cj;_1, and J3;=Cy_;. If the three 
parts are now considered as a whole, the corre- 
sponding polynomial, which represents the to- 
tality of structures with a bond between orbitals 
1 and 2), is the product J;/2J/3. This follows from 
the fact that there can be no bonds between any 
two of the three parts, since such bonds would 
necessarily lead to noncanonical structures. The 
law for combining such noninteracting systems 
is then formally identical with that for multi- 
plying polynomials. We now consider the canon- 
ical structures in which the orbital 1 is bonded 
with some further definite orbital, say 2k. As in 
the previous case, the associated polynomial is 
equal to K,K2K3, with K,=1 or 2, Ke=C;,_1, 
and K;=Cy_,. Obviously none of the structures 
of this second group will occur also in the first 
group, and consequently the totality of structures 
in which the orbital 1 is bonded to either 27 or 2k 
will be represented by the sum J\J2J3+K1K2K3. 
This same procedure can now be repeated until 
all of the possible canonical structures have been 
taken into consideration. The resulting poly- 
nomial is then 


N 
Cy =VSiJ2J3, 


j=1 


which is just the recursion formula given 
originally. 
Since we have adopted the convention that 
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Cy=1, and since we see at once that C,;=1, we 
can find the remaining polynomials in order. 
The first few members of the series are: 


C,=1, 

Co=1+ Z, 

C3=1+ 324+ 2, 

Cy=1+ 624+ 6224+ ~ 2, 
C;=1+102+ 2027+ 102°+21. 


In general : 


N-1 


N\N-1)! 
Cr= 7%, -— ——s 
i=0 7 1(j+1) (N-—7) (N-j-1)! 


i 





1 ‘-1/N N 
EAs as)” 
N i-o\ G7 \j+i 


II. THe MoLecuLe Consists oF A RING 


For a ring of 2N orbitals, we can carry through 
the same sort of procedure as in the case of the 
chain. The result is then 


N-1 
Ry =2Cy_-1+2 p> Ci-1C nw} 


j=2 
This can be expressed more conveniently as 
Ry =Cy+(1—2)Cw-1, 


in which form it is seen to be a special case of 
the general rule given below in III. The first few 
members of the series are: 


R2=2, 

R3=2+ 32, 

R,=2+ 82+ 42’, 
R;=2+152+ 2027+ 52°. 


III]. THE MOLECULE CONTAINS SEVERAL RINGS 


It is convenient to represent the polynomial 
associated with a given molecule by means of a 
geometrical figure, which is constructed in the 
following way : After the orbitals are arranged in 
the circle, or polygon, lines are drawn connecting 
each pair of orbitals which are adjacent in the 
actual molecule, that is, each pair of orbitals for 
which the corresponding exchange integral is to 
be set equal to a. Thus we have, for example: 





G. W. 


[| 


() 
\“F 


and so on. 

Now let us consider the relation between two 
figures which differ only in that one contains a 
single line which is absent from the other, as for 


| ta di hth 
example: , corresponding to naphtha- 
Vr a 


PSS 


lene, and | |, corresponding to cyclo- 
Ws 

decapentane. For such cases, there is a general 

rule: If a figure, A, differs from a second, B, by 

having a single additional line, and if this line, 

together with its end points, divides A into two 

parts, C and D, then A=B+(1i-—z)CD. The 
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derivation of this equation follows from the same 
considerations as were used in I. Thus we con- 
sider the totality of the structures of each 
molecule in which there is a bond along this 
added line. In A these are represented by the 
product CD, and in B by the product 2CD. 
Structures which do not contain this bond are 
not affected in degree of excitation by the added 
line. This completes the proof. The theorem, 
when applied to naphthalene, gives 


oY~N se J\(N\N) 
| = + (1—2);| IT 
f~/ ANS NY] 


=R;+(1-—2)C?? 
=3+162+ 1927+ 42°. 


We find, then, that naphthalene has 3 unexcited, 
16 first excited, 19 second excited, and 4 third 
excited structures, in agreement with the direct 
evaluation by Pauling and Wheland.! 


A slightly more complicated example is afforded by biphenyl.°® 


2 J 


LH OO Hal | 


~ 


a a," 
= i 


+0-){C 


— Cet+(1 —2z)C3Co+(1 —2)CoR3 
= 4+ 2124+ 482?4+ 4423+ 142!+ 2°. 


This result has been checked by writing down all 132 canonical structures, and counting the number 


of each degree of excitation. 
A few further examples are given below: 


Phenylethylene 


Sym. diphenylethylene 


Unsym. diphenylethylene 


5 Here, and in the following, we do not arrange the 
orbitals actually in a circle or polygon, but in whatever 
type of figure is most convenient for the purpose at hand. 
In cases in which there may be some question, the orbitals 


=2+62+52?+32', 


=4+ 332+ 1132+ 1652°+-942!+ 1925+ 28, 


= 44+ 2424+-952?+ 1322°+ 1342!+ 3325+ 72°, 


will be numbered, to indicate the sequence which they 
would follow if they had been arranged in the polygon. 
It is absolutely essential that the same sequence be 
maintained throughout a given calculation. 
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Anthracene = 4+ 482+ 1502?+ 1632+ 5824+ 625, 


Phenanthrene =5+472+ 1482?+ 1652*+ 5924+ 525, 


9 ,. 
Naphthacene =5+1102-+6492?+ 1556s°+ 164624+750z5 
' +1382°+82". 


1,2 Benzanthracene =7+1122+6422?+ 15512°+ 1654214 7545 
+1352°+-72’, 


3,4 Benzophenanthrene = 8+ 1162+ 6362? + 15422°+ 166324+ 76025 
+1312°+62", 


Chrysene =8+11724+ 6342?+ 15412°+ 166724!+ 75925 
+1292°+-72’, 


Triphenylene =9+1172+ 6452°+ 155723 + 165924+ 74425 
' +1252°+ 6:2’, 


Terphenyl =8+812+ 3842?+ 10002?+ 152824+- 127925 
f +5032°+ 7627 +328. 


In the resolution of several of the above figures, it is convenient to make use of a further rule which 
can easily be derived from the one above: Let D>», , be the polynomial corresponding to two straight 
chains, of and q orbitals, respectively, which are not connected to each other. Then 


Dy, (= Cort get (2—- 1)Dp-1, q-1 


(p+4q—|p—aq|)/2 


= } (2—1)*C(p4 /2-%- 


k=0 
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Ma = 122+ 3227+ 622°+ 202'+ 62°. 


IV. THe MoLecuLE CONTAINS AN UNPAIRED 
ELECTRON 


This type of problem arises in the treatment of 
free radicals. The procedure is simply to intro- 
duce a “phantom orbital” X, containing a 
“phantom electron,” which is formally bonded 
to the unpaired electron of the molecule.'! The 
exchange integrals between this phantom orbital 
and the real orbitals are set equal to zero. As 
before, it is desirable first to obtain expressions 
for straight chains, and we find at once that the 
polynomial Fy, corresponding to the chain of 
2N-—1 orbitals, is given by the expression 


Fy = (1/2)Den-1, 1= (1/z)[Cv+ (g— 1)Cy_1 ]. 


(The factor 1/z occurs since the bond to the 
phantom orbital X, although ineffective, is not 
considered to increase the degree of excitation.) 





A few other results are given below: 


a-naphthylmethyl 


6-naphthylmethyl 


Biphenylmethyl 


APPLICATION TO THE CALCULATION OF 
MOLECULAR ENERGIES 


The author has shown recently® that for a 
hydrocarbon molecule which does not contain 
any odd membered rings the energy which results 
from the simplified Heitler-London procedure of 

6G. W. Wheland, J. Chem. Phys. 3, 230 (1935). 


= Cet (2—1)C5+(2—1)?C.+ (sz —1)°C34+- (2-— 1)*C,+(2—1)°Ci, 





The first few members of the series are: 


F,=1, 

F.=2, 

F3=3+ 2z, 

Fy=4+ 82+ 22?, 
F,;=5+202+ 152?+22'. 


It is interesting to observe that the coefficients 
k; appearing in Fy are the same as those appear- 
ing in Ry, but that their order is reversed. 

In the treatment of molecules containing rings, 
the same procedure can be used as before. For 
the phenylmethy! radical, for example, we have 


2 3 
« 
4 


x 


a _vf7 Vex 
Cte ft * 
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=F,4+(1 —2)C,F,=5+82+2". 


= 10+462+582?+ 182%, 


=9+502+ 5422+ 1823 +24, 


= 13+ 662+ 1482? + 1472°+ 5124+ 42°, 





Pauling and his co-workers! can be expressed as 


- DL RAR RNR 
sia DLrdegr ~’ 


where ge is the number of structures of type X 
which appear with the coefficient az in the com- 
plete eigenfunction for the molecule, and mz '5 
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CANONICAL STRUCTURES FOR HYDROCARBONS 


the number of effective bonds in any one of the 
structures R. Q and a are, as usual, the Coulomb 
integral and the single exchange integral, re- 
spectively. If all the structures of the ith degree 
of excitation are assumed to have an average 
coefficient a;, then the above equation becomes 


> iaiki(b—1) 
Lak ™ 


where b is equal to the number of effective bonds 
in the unexcited structures. (If there are no 
phantom orbitals, b= N=one-half the number of 
orbitals.) 

If some method were available for estimating 
the coefficients a;, the calculation of the energy 
of any molecule of the present type would be 
quite simple. The rule given in the previous 
paper,® that a;=(0.3)‘ is fairly satisfactory for 
small molecules, but it leads to too small values of 
the energy (in absolute magnitude) when NV 
becomes greater than 5. It also fails to allow for 
variations of the coefficients between molecules 
of the same size, but of quite different structure. 
A rather better, but equally empirical rule can 
be put forward tentatively : 


a;=(0.04—0.93/N+1.05B/N?)i, 


where B is equal to the number of pairs of 
orbitals for which the corresponding exchange 
integral is a. Some results obtained by this 
method are given in Table I. It will be seen that 
they are really quite satisfactory. It is doubtful, 
however, if this procedure could be safely applied 
to molecules of widely different type from those 
considered here, since there is no assurance that 
the same rule for estimating coefficients would be 
at all satisfactory.’ In particular, the assumption 


W=0+ 


7One readily sees, in fact, that the rule is no longer 
very good for the free radicals. 


TABLE [, 








Energy 
Approx. Correct* 


Q+2.6la 
3.3la 
4.05a 
3.30 
5.66a 
5.67a 
4.93a 
7.36 
4.08a 
5.59a 
5.65a 
7.15a@ 
7.21la@ 
7.24c 
7.240 
7.262 
1.74a 
2.45a 
3.21a 


Compound a 


0.43 
0.33 
0.27 
0.33 
0.23 
0.23 
0.26 
0.20 
0.32 
0.25 
0.25 
0.21 
0.21 
0.21 
0.21 
0.21 
0.36 
0.31 
0.27 





0+2.6la 
3.30a 
4.03a 
3.31ae 
5.59ae 
§.55a° 
4.87a¢ 
7.0908 
4.04a 
(5.53a)> 
(5.69a)» 


Benzene 
Cyclooctatetrane 
Cyclodecapentane 
Phenylethylene 

Sym. Diphenylethylene 
Unsym. Diphenylethylene 
Biphenyl 

Terphenyl 
Naphthalene 
Anthracene 
Phenanthrene 
Naphthacene 

1,2 Benzanthracene 
3,4 Benzophenanthrene 
Chrysene 
Triphenylene 
Butadiene 

Hexatriene 
Octatetrane 


(7.35) 


1.73a 
2.48a 
3.22a¢ 








8L. Pauling and J. Sherman, J. Chem. Phys. 1, 679 (1933). G. W. 
Wheland, J. Chem. Phys. 3, 230 (1935). 

b These values were calculated from the observed resonance energies 
(L. Pauling and J. Sherman, J. Chem. Phys. 1, 606 (1933)) with a= --1.5 
e.v. In view of the really surprising agreement which has been found, in 
all other cases, between the theoretical and the observed resonance 
energies (cf. ref. 1), this seems a safe method for estimating the order of 
magnitude of the former. 

¢In these calculations, Pauling and Sherman used an approximate 
method of solving the secular equation. The procedure (for a justifica- 
tion of which see J. Sherman, J. Chem. Phys. 2, 488 (1934)) consisted in 
neglecting the second and higher excited structures and in arbitrarily 
equating coefficients. As a result, the coefficients of a appearing here are 
somewhat too low. 


that a; for a given molecule can be represented 
by the ith power of a constant quantity has 
little if any theoretical justification. The esti- 
mated coefficients may, therefore be very seri- 
ously in error, and the energy may come out 
approximately correct only in consequence of a 
cancelation of errors. Be this as it may, however, 
the method seems to be quite a satisfactory one 
from an empirical point of view, and can be 
applied with some assurance to the treatment of 
molecules which do not differ very greatly, either 
in size or in type, from those considered in 


Table I. 
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The Thermal Conductivity of Polyatomic Gases. Erratum 


A. R. UBBELOHDE, Department of Thermodynamics, Clarendon Laboratory, Oxford 
(J. Chem. Phys. 3, 219 (1935)) 


AGE 222, 6th line from the bottom: Should read, ‘‘a maximum addition of 
0.02 to F”’’ in lieu of ‘ta maximum of 0.02 to F’’. 





Molecular Arrangement and X-Ray Diffraction in Ionic Solutions. Errata 


J. A. Prins, Natuurkundig Laboratorium der Rijks- Universiteit, Groningen, Holland 
(J. Chem. Phys. 3, 72 (1935)) 


Y article on this subject having been 
printed without author’s proof, the follow- 
ing inaccuracies occur: 


Page 77, end second paragraph §7. To be added: 
It should be remarked that at a larger 
distance from the ions a different and more 
subtle rearrangement of the water mole- 
cules may take place. It is this rearrange- 
ment (‘‘water II—water I’’) which has 
been considered by Bernal and Fowler in 
connection with the term ‘‘structural tem- 
perature.” It will presumably be of im- 
portance in moderately diluted solutions. 
In our case, however, of strong solutions, 
where every water molecule is rather close 


to an ion, I think the rearrangement con- 
sidered above (‘‘water II—water III’’) is 
more important. 

Page 78, paragraph headed by “If the case of a 
1: 10 Lil solution.” To be read: class III 
mentioned in §6, instead of class II, men- 
tioned in §4. 

Page 79, beginning. To be read: class II ($5) 
instead of class III. 

Page 80. To be read: 1.8 instead of 2.8A. 

Page 80. To be read: (a+l1+6+1)/2 instead of 
a+l+b+. 

I use this occasion to remark that a more 
quantitative study of some points considered in 
this article will be published in Physica (Den 
Haag) in collaboration with Dr. R. Fonteyne. 
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This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section must reach the office of the Managing Editor 


The Raman Spectrum of Orthophosphoric Acid 


The Raman spectrum of orthophosphoric acid has been 
investigated by Nisi.! He found frequency differences of 
356, 498 and 911 cm for 84 percent (by weight?); 512 and 
906 for 50 percent; and 505 and 898 for 25 percent. The 
writers have made a study of H;PO, at concentrations 
ranging from 50 mole percent (the highest available) to 10 
percent. The results agree with those of Nisi. 

One strong line appears at all concentrations. Measure- 
ments of it are as follows (with the excitation by 4358A 
accurate to about 2 cm™’): 


Mole percent Excited by 4358 Excited by 4047 


50 22,029 909 
23,798 907 


40 22,033 905 
23,801 


30 22,037 
23,807 


22,042 
23,807 


22,043 895 
23,812 893 


When the product of exposure time and mole percent is 
made constant, this line appears to have the same intensity 
at all concentrations. A shift in frequency of about 14 cm=! 
is clearly indicated. 

Other lines which appear at 50 percent are: 22,448 (490), 
22,579 (359), 23,607 (909), and 24,212 (493). 490 and 359 
cm™ are broad and weak. The 23,607 line is probably ex- 
cited by 4078A. These weak lines could not be measured at 
lower concentrations, except 23,624 (892) at 10 percent. 

Myron A. JEPPESEN 
RAYMOND M. BELL 
Department of Physics, 
Pennsylvania State College, 
April 22, 1935. 


!Nisi, Jap. J. Phys. 5, 119 (1929). 


The Vitreous State 


Inarecent article,! Zachariasen has commented my paper 
on the vitreous state, published in the January number of 
this Journal.? In the following, the comments of Zachariasen 
are briefly discussed. 


not later than the 15th of the month preceding that of the 
issue in which the letter is to appear. No proof will be sent 
to the authors. The usual publication charge ($3.00 per 
page) will not be made and no reprints will be furnished free. 


My statement that Zachariasen’s viewpoint is only 
adaptable to cases in which the oxygen polyhedra form 
three-dimensional networks, is based on the fact, that 
Zachariasen only treats this case. The requirement of a 
three-dimensional network is the fourth of his rules which 
he thinks must be fulfilled if an oxide is to form a glass. 
Later, this rule is also generally applied to oxides, without 
considering the possibility that the corresponding crystal 
structure could be characterized by a network of oxygen 
polyhedra, extended in less than three diinensions. Also in 
his summary, Zachariasen speaks only of three-dimensional 
networks. 

Zachariasen further says that it is “no real progress to 
present a theory for the structure of glass in terms of the 
equally unknown structure of the liquids.” I think, in con- 
trary, that this is the most logical way of attacking the glass 
problem. We know from several analogies that the liquid 
and vitreous states must be closely connected and we then 
make use of conceptions of the liquid state which are so 
general and obvious that they are very likely to be true. I 
do not think this is more dangerous than direct compari- 
sons between the vitreous and the crystalline state. 

The fact that I have not given any detailed description 
of the structure of the fragments, possible in a glass, is 
caused by the multitude of structures existing in different 
glasses and by the variation which can be expected for the 
different fragments in one definite glass. I have pointed out 
the importance of definite coordinations within the frag- 
ments and I think a consideration of this fact is sufficient 
for making clear the general structure of a fragment. 

It is a matter of course, that a glass, built according to 
my theory, must be optically clear as the constituting 
fragments or groups are small compared to the wavelength 
of light. 

It is my opinion that it is the tendency of certain atoms 
to attain a certain oxygen coordination and the strength 
with which oxygen is bound to these atoms which is the 
main cause of the glass forming tendency of vitreous oxides. 
As Zachariasen remarks, I added in the summary that “the 
glass forming tendency . . . seems to be especially large 
when the corresponding crystals contain radicals or mole- 
cules unlimited in one or more dimensions” but I do not 
admit that this statement is too general when it is read in 
connection with the preceding sentence which implies that 
these groups must be contained in the melt and kept to- 
gether with strong (probably in many cases homopolar) 
forces. In this way I explained why SiO, can form glass and 
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BeO not. The fact that TiO, is not known in the vitreous 
state is also explained by the more basic character of this 
oxide, which is connected with weaker bonds between Ti 
and O. 

Concerning H;PO,, I explained its glass forming tend- 
ency by assuming larger groups than (PO,)~* to exist in the 
melt but I did not discuss the question whether these groups 
are formed by polymerization of (PO,)~* radicals only or if 
hydrogen was necessary at the construction. 

I consider the glasses as supercooled liquids in which the 
supercooling tendency is caused by groups in the melt 
which can be added with difficulty directly to the lattice 
and within which the atoms are kept together by so strong 
forces as to make their separation difficult. My theory of 
the structure of glasses is based on this assumption and the 
problem of the glass forming tendency of a substance is 
brought back to the question whether such groups are to be 
expected in the melt. Zachariasen, on the contrary, has 
assumed a certain type of structure, characterizing the 
vitreous state. Starting from the corresponding crystal 
structures, he then discusses the possibility of constructing 
the vitreous structure. 

GuNNAR HAGG 

Institute of General and Inorganic Chemistry, 

University of Stockholm, 
April 2, 1935. 


1 W. H. Zachariasen, J. Chem. Phys. 3, 162 (1935). 
2G. Hagg, J. Chem. Phys. 3, 42 (1935). 


The Magnetic Properties of the Phenanthroline Ferric 
Complexes 


We have investigated the magnetic susceptibilities of 
the complex ions formed by the base ortho-phenanthroline 
with ferric ion. Blau, who discovered some forty years ago 
the extraordinary stability of the complex metallic ions 
formed by this base, reported that two ferric complexes are 
formed, which are not directly jnterconvertible. One of 
these is intensely blue in color; can be obtained only 
by oxidation of the ferrous complex; and has, as the 
analysis of several solid salts demonstrate, the formula 
Fe(CizHsN2)3***. The other is brown, and is formed by 
the direct reaction of ferric salts with phenanthroline. 
We have obtained a solid brown salt whose analysis cor- 
responds to the formula [Fe(Ci2HsNe)2(H2O)(OH) JCle, 
but other compounds are undoubtedly also present in the 
solutions obtained by addition of phenanthroline to solu- 
tions of ferric salts. As is indeed to be expected of complexes 
as stable as these, their formation from ferric ion is ac- 
companied by a very considerable decrease in magnetic 
susceptibility. According to the theories of Lewis and 
Pauling the linkages concerned are therefore to be consid- 
ered as of the shared electron rather than of the ion-dipole 
type. From measurements in aqueous solution the atomic 
susceptibility of the iron atom in the blue complex is 
2446 X10-* which is practically identical with the value 
observed with ferricyanide ion and generally accepted as 
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the normal one for coordinatively bound ferric ion. The 
corresponding value for the brown complex is, however, 
only 826X10-*. The low value for the magnetic moment, 
1.4 Bohr magnetons, to which this leads is perhaps to be 
attributed to the existence in a single molecular grouping 
of two or more iron atoms so oriented that their moments 
partially neutralize each other. There are indeed other 
indications, now under investigation, that the brown com- 
plex has a higher molecular weight than would correspond 
to a single iron atom per molecule. 
GEorRGE H. WALDEN 
Louis P. HAMMETT 
ALLISON GAINES, JR. 
Department of Chemistry, 
Columbia University, 
May 6, 1935. 


The Decomposition of Hydrogen Peroxide by the Irradi- 
ation of Its Aqueous Solution with X-Rays* 


Gas-free (unbuffered) solutions of hydrogen peroxide 
were irradiated with x-rays and the resulting decomposition 
determined by electrometric titration! and by gas analysis. 
None or a very small amount of hydrogen is produced, the 
decomposition proceeding essentially according to: 2H,O» 
=2H,0+02. Over the range of concentrations (0.0001 to 
0.1 moles H,0,/1000 cc) and x-ray intensities (3 to 15 
r/sec.) used, the number of moles of hydrogen peroxide 


decomposed per unit of dosage increases as the square root 
of the concentration and decreases as the square root of the 
x-ray intensity. (Fig. 1.) For a special set of conditions, the 
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Fic. 1. Decomposition of hydrogen peroxide as a function of: (A) 

Square root of concentration (C in micromoles per 1000 cc) for an 


intensity of 9.5 r/sec. (B) Inverse square root of intensity (/ in r/sec.) 
for a concentration of 1.0 millimole H2O2 per 1000 cc. 


dependence of the decomposition** on the temperature was 
determined between 0° and 50°C. (Fig. 2.) The temperature 
coefficient shows a remarkable rise with the temperature; 
its value, between 15°C and 25°C is 1.18.7 

A comparison with results?-* obtained for the decomposi- 
tion of hydrogen peroxide with light is of interest. The 
evidence is not wholly clear, but the increase of the 
quantum efficiency with the square root of the concentra- 
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Fic. 2. Decomposition of hydrogen peroxide as a function of the 
temperature for an x-ray intensity of 3.0 r/sec. and a concentration of 
1.0 millimole H2O2 per 1000 cc. 


tion (except when the concentration is very large)” * and 
with the inverse square root of the intensity, appears con- 
sistent with the experimental data. This would be in agree- 
ment with the present results if we assume that the de- 
composition in the case of x-rays is not due to the direct 
activation of the hydrogen peroxide molecules by the x-ray 
electrons, but that the primary process is the activation of 
water molecules! 3 followed by the transfer of energy 
from these activated molecules to the hydrogen peroxide.tt 
On this assumption, the quantum efficiency (in the case of 
light) should be compared to the number of moles of 
hydrogen peroxide decomposed per unit of x-ray dosage. It 
should also be noted that if the decomposition were due to 
the direct activation of the hydrogen peroxide molecules by 
the x-ray electrons, in the extreme case as many as 10,000 
molecules of hydrogen peroxide would be decomposed for 
each hydrogen peroxide molecule activated, an unlikely 
number in view of the relatively short chains observed in 
the case of light. 
HuGo FRICKE 
Walter B. James Laboratory for Biophysics, 
The Biological Laboratory, 
Cold Spring Harbor, N. Y., 
March 30, 1935. 


*1I am indebted to Dr. E. R. Brownscombe for assisting with these 
measurements. . 

** The decomposition of nonirradiated blanks was corrected for, and 
amounted, at 50°C, toa few percent of that produced during irradiation. 

1 For the decomposition by light, Allmand and Style, reference 2, 
found the temperature coefficient (at room temperature) to be: 1.42 
for A=3650A, and 1.28 for \=2750A. 

tt On this assumption, in our experiments the number of molecules 
of hydrogen peroxide activated would be expected!, 9-13 to be of the 
order of 10'3/cc/sec. With light, intensities producing from 103 to 1016 
activations/cc/sec. have been used. 

1 Fricke, J. Chem. Phys. 2, 556 (1934). 

*Allmand and Style, J. Chem. Soc. (1930), 596, 606. 
one Comptes rendus 151, 1040 (1910); 156, 1601, 1758, 1879 
wi and Wurmser, Comptes rendus 156, 1012 (1913); 157, 126 


* Kornfeld, Zeits. f. wiss. Phot. 21, 66 (1921). 

‘Winther, Trans. Faraday Soc. 21, 459 (1925). 

: Heidt, J. Am. Chem. Soc. 54, 2840 (1932). 

; Qureshi and Rahman, Zeits. f. physik. Chemie 36, 664 (1932). 
Fricke and Morse, Phil. Mag. 7, 129 (1929). 

" Fricke and Brownscombe, J. Am. Chem. Soc. 55, 2358 (1933). 
- Fricke and Brownscombe, Phys. Rev. 44, 240 (1933). 

Fricke and Hart, J. Chem. Phys. 2, 824 (1934). 

* Fricke and Hart, J. Chem. Phys. 3, 60 (1935). 
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The Oxidation of the Nitrite to the Nitrate Ion by the 
Irradiation of Its Aqueous Solutions with X-Rays 


Gas-free solutions containing the nitrite ion were ir- 
radiated with x-rays and were analyzed for hydrogen and 
oxygen; and for nitrite by electrometric titration with 
potassium permanganate. The nitrite ion is oxidized to the 
nitrate ion and hydrogen liberated in the equivalent 
amount. The amount of nitrite oxidized is independent of 
the concentration of the nitrite ion from 0.05 to 100.0 milli- 
moles/1000 cc and of the hydrogen ion concentration from 
pH =2.0 to 11.0 and equals 0.55 micromole/1000 r/1000 
cc. It is concluded that the transformation is due to water 
molecules activated by the x-rays, according to: 

NO2- + (H20) act = NOs~ + He. 
HuaGo FRIcKE 
Epwin J. Hart 
Walter B. James Laboratory for Biophysics, 
The Biological Laboratory, 
Cold Spring Harbor, Long Island, N. Y., 
May 13, 1935. 


The Far Ultraviolet Absorption Spectra and Ionization 
Potentials of the Ethyl Halides 


The absorption spectrum of ethyl iodide is essentially 
similar to that of methyl iodide! and starts at about the 
same wavelength. Two electronic series containing a large 
number of terms appear. The upper members can be ex- 
trapolated in the form of Rydberg series to ionization 
potentials of 9.295+0.005 and 9.885 +0.005 volts. Thus a 
difference of about 0.6 volt is maintained in going from 
methyl to ethyl iodide though the ionization potentials of 
the latter are diminished by 0.2 volt relative to the cor- 
responding ones in methyl iodide. 

The first resonance bands of ethyl bromide start around 
1770A and thus suffer a short wavelength shift of 15A 
relative to the corresponding bands in methyl bromide. 
They are much more diffuse but two series limits can still 
be observed. These can be extrapolated to the ionization 
potentials 10.24+0.01 and 10.56+0.01 volts corresponding 
to a diminution of 0.25 volt relative to the values obtained 
for methyl bromide.’ Part of the difficulty of observing 
these limits is a background of continuous absorption 
arising from excitation in the CC bond. 

In the case of ethyl chloride only the first few bands are 
sharp. These are shifted to long wavelengths by about 
0.36 volt relative to the analogous bands in methy] chloride. 
The electronic difference of 0.08 volt is maintained and it is 
estimated that the ionization potentials of ethyl chloride 
are about 0.38 volt less than those given for methyl 
chloride. 

The above data provide considerable support for the 
suggestion put forward by the author® that the dipoles 
of methyl and ethyl groups have a strong influence in 
reducing the ionization potentials of electrons in adjacent 
bonds. R. S. Mulliken‘ has shown that the spectra of the 
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alkyl halides are due to the excitation of nonbonding np 
electrons from the halogen atom. It is here assumed that 
the additional strength of the C2H; dipole over the CH; 
dipole is responsible for the diminution of the ionization 
potentials of the ethyl relative to the methyl halides. 
The increasing effect as we go from the iodide to the chloride 
is explained by the decreasing combination radii of the 
halogen atoms. This brings the dipole closer to the electron 
considered and thus increases its effect. 

This dipole hypothesis has also been useful in explaining 
the deviations of the first few electronic states from Ryd- 
berg formulas which occur if the spectrum starts at fairly 
long wavelength (i.e., around 2000A). These deviations 
are negligible in the case of C2H2, H2S and CH;Cl which 
start around 1500-1600A. In these cases the radii of the 
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orbits of the first excited states are so large that the effect 
of the dipole disappears after the first jump. In methy] 
and ethyl iodide the first and second excited states are still 
strongly affected by the dipole. The hypothesis can also be 
applied to give a satisfactory explanation of the magnitudes 
of the long wavelength shifts of the first resonance bands 
of the ethyl halides which appear to be a little anomolous 
at first sight. 
W. C. Price 
The Rowland Physical Laboratory, 
The Johns Hopkins University, 
May 14, 1935. 
1W. C. Price, Phys. Rev. 47, 419 (1935). 
2W. C. Price, Phys. Rev. 47, 510 (1935). 


3 W. C. Price, J. Chem. Phys. 3, 256 (1935). 
4R. S. Mulliken, Phys. Rev. 47, 413 (1935). 





